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ABSTRACT 

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini 
OBI molecular cloud complex, and targeted NH3 observations of the BGPS sources. When paired with 
molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such 
as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS 
sources above 5a = 0.37 Jy beam~^ with corresponding 5a detections in the NH3(1,1) transition. 
Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 
1.12 pc, and a mean gas kinetic temperature of 20 K for the sample of 34 BGPS sources with detections 
in the NH3(1,1) line. The observed NH3 line widths are dominated by non-thermal motions, typically 
found to be a few times the thermal sound speed expected for the derived kinetic temperature. We 
calculate the mass for each source from the millimeter flux assuming the sources are isothermal and 
find a mean isothermal mass within a 120" aperture of 230 ± 180 Mq. We find a total mass of 8,400 
M0 for all BGPS sources in the Gemini OBI molecular cloud, representing 6.5% of the cloud mass. By 
comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source 
size calculated from the NH3 line widths, we find a mean virial parameter {M^r / Miso) of 1.0±0.9 for 
the sample. We find mean values for the distributions of column densities of 1.0 x 10^^ cm^^ for H2, 
and 3.0 x lO^'' cm^^ for NH3, giving a mean NH3 abundance of 3.0 x 10^* relative to H2. We find 
volume-averaged densities on the order of 10'^ — 10^ cm^^. The sizes and densities suggest that in 
the Gem OB 1 region the BGPS is detecting the clumps from which stellar clusters form, rather than 
smaller, higher density cores where single stars or small multiple systems form. 

Subject headings: stars: formation — ISM: individualGem OBI — ISM: dust — ISM: clouds — radio 
lines:ISM 



1. INTRODUCTION 

The availability of millimeter bolometer arrays has 
made possible large-scale, blind surveys of the cold dust 
most closely associated with star formation (e.g., Enoch 
et al. 2007; Motte et al. 2007; Aguirre et al. 2010; 
Schuller et al. 2009). Such blind millimeter continuum 
surveys can identify the clumps where massive stars and 
clusters are born without the requirement of a signpost 
of massive star formation, such as maser emission (e.g., 
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Palla et al. 1991; Plume et al. 1992; Plume et al. 1997; 
Beuther et al 2002; Mueller et al. 2002; Shirley et al. 
2003), radio continuum emission (e.g., Ramesh & Sridha- 
ran 1997; Sridharan 2002), infrared sources (e.g., Moli- 
nari et al 2000; Kumar et al. 2006; Kumar & Grave 
2007; Robitaille et al 2008), or infrared colors (Wood & 
Churchwell 1989). Infrared dark clouds (IRDCs; regions 
of high density seen in absorption against the diffuse mid- 
IR background which are thought to represent the earli- 
est stages of massive star formation) provide another way 
to find regions of massive star formation (Rathborne et 
al. 2006; Jackson et al. 2008), but only if they lie in front 
of extended mid-IR emission. However, IRDCs are seen 
in emission at mm wavelengths, allowing the BGPS to de- 
tect the earliest stages of massive star formation at larger 
distances where a significant mid-IR background is ab- 
sent or too much foreground mid-IR emission is present. 
Galaxy-wide millimeter continuum surveys, in con- 
junction with distance information, will allow us to deter- 
mine the physical properties of a vast number of sources 
and begin answering some important questions. What 
is the impact of environment on star formation? How 
do the physical properties of clumps, such as size, mass, 
and mean density vary as a function of Galactocentric 
radius and proximity to a spiral arm? Answering these 
questions will provide a greater understanding of star for- 
mation in our own Galaxy. Also, by characterizing the 



distribution of dense gas clumps in our own galaxy, we 
can better understand the aggregate emission of clumps 
in distant galaxies where the dense gas structures cannot 
be resolved. 

In this paper, we refer to clumps as distinct bound re- 
gions within molecular clouds that will form entire stellar 
clusters, and cores as the regions within clumps that will 
form a single star or small multiple system (e.g., Williams 
et al. 2000). Clumps can be characterized with masses 
ranging from 50 — 500 Mq, sizes ranging from 0.3 — 3 
pc, mean densities ranging from 10'^ — 10^ cm^^, and gas 
temperatures ranging from 10 — 20 K, while starless cores 
can be characterized by masses from 0.5 — 5 Mq, sizes 
of 0.03 — 0.2 pc, mean densities of 10** — 10^ cm~^, and 
gas temperatures between 8 and 12 K (Bergin & Tafalla 
2007). 

The Bolocam Galactic Plane Survey (BGPS) is one 
of the first millimeter surveys of the northern Galactic 
plane (Aguirre et al. 2009), and has cataloged 8,358 
sources (Rosolowsky et al. 2009). The BGPS consists 
of two parts: a blind survey of the inner Galaxy, and 
a targeted survey of well-known star-forming regions in 
the outer Galaxy. The inner Galaxy survey coverage is 
continuous over —10° < I < 90.5°, and the outer Galaxy 
coverage includes IC1396, the W3/4/5 region, 4 square 
degrees towards the Perseus Arm (near NGC 7538), and 
the Gemini OBI molecular cloud. With a large range in 
Galactic longitude, the BGPS provides the opportunity 
to study star formation in drastically different environ- 
ments. 

Millimeter continuum emission detected in the BGPS 
can arise from all objects along a line of sight, from a few 
hundred pc to the other side of the Galaxy; distances 
to BGPS sources are not known a priori. With recent 
improvements to kinematic models of the Galaxy (e.g., 
Pohl et al. 2008) and VLBI parallax distances (Reid et 
al. 2009), radial velocities can yield kinematic distances. 
Molecular spectral lines also provide information regard- 
ing the kinematics within each source. Ammonia (NH3), 
in particular, is an excellent tracer of dense gas which also 
provides information regarding the line of sight velocity 
dispersion, kinetic gas temperature, column density and 
mass (e.g.. Ho & Townes 1983; Zinchenko et al. 1997; 
Swift et al. 2005; Rosolowsky et al. 2008; Foster et al. 
2009). Pairing the 1.1 mm continuum data with NH3 ob- 
servations allows complete characterization of both the 
gas and dust properties (e.g., Rosolowsky et al. 2008; 
Foster et al. 2009). 

In this paper, we present the results of an NH3 survey 
of the BGPS sources in the Gemini OBI molecular cloud. 
In ^ we provide an overview of the structure and pre- 
vious observations of the Gemini OBI molecular cloud. 
We briefly describe the BGPS and NH3 surveys and data 
reduction in fj3l and present the parameter estimation 
methods in 21 In El we discuss the basic results of the 
BGPS and NII3 in the Gemini OBI molecular cloud, and 
in fj6]we discuss the derived physical properties in detail. 
We compare our results to previous studies in §7.1) and 
explore the effects of distance on the derived properties 
in ^ 37. 21 Finally, we present a summary in f}8l In fu- 
ture work, we will present GBT-t-BGPS observations of 
sources at low Galactic longitudes. 

Studying the Gem OBI region is particularly useful be- 
cause the BGPS sources in this region are located in the 



same molecular complex and thus have similar distances 
(Carpenter et al. 1995a; Carpenter et al. 1995b). We are 
able to investigate the physical properties of this subset 
of BGPS sources without the ambiguity introduced by a 
large range of distances seen for the sample as a whole. 
This paper will serve as a Galactic anti-center compar- 
ison for a study of NII3 observations of BGPS sources 
at a range of Galactic longitudes in the inner Galaxy 
(M. K. Dunham et al., in prep.). 

2. GEMINI OBI MOLECULAR CLOUD 

Gemini OBI is an OB association centered near {I, h) = 
(189°, 1°) which was first identified by Morgan, Whit- 
ford, & Code (1953). It became known as Gem OBI 
when Humphreys (1978) compiled a list of the brightest 
OB stars in the Milky Way, cataloging 20 OB stars as 
belonging to the association. Many groups have stud- 
ied the association and found distances ranging from 1.4 
kpc to 1.9 kpc (Crawford et al. 1955; Hardie et al. 1960; 
Georgelin et al. 1973; Humphreys 1978). There is also a 
significant spread in age amongst the association mem- 
bers, with ages ranging from 2 Myr (Grasdalen & Car- 
rasco 1975) to 10 Myr (Barbaro et al. 1969). 

The Gemini OBI molecular cloud was first identified 
as an enhancement of emission in the Columbia ^^CO 
survey of the Galactic plane (Dame et al. 1987), and 
spans a 6° X 6° area on the sky (Huang & Thaddeus 
1986). There are many star-forming regions in the Gem 
OBI molecular cloud, including well-known HII regions 
S247, S252, the group S254, S255, S256, S257, and S258 
(hereafter referred to as S254-258), S259, S261 (Sharpless 
1959), and BFS 52 as well as a supernova remnant, IC 
443. Reid et al. (2009) used the VLBA to measure the 
parallax of methanol masers in S252, and found a dis- 
tance of 2.10 ± 0.026 kpc, which we adopt in this study. 

Studies of ^^CO and ^^CO in the Gem OBI molec- 
ular cloud show abundant arc and ring-like filamentary 
structure (Carpenter et al. 1995a; Kompe et al. 1989). 
Carpenter et al. (1995a) mapped 32 square degrees of 
the Gem OBI molecular cloud in ^^CO(J = 1 — 0) and 
13C0( J = 1-0) with 50" sampling and 45" and 50" res- 
olution respectively. They suggest that the filamentary 
structure is comprised of molecular material which has 
been swept-up by the expanding HII regions and wind 
blown bubbles. They conclude that the high column 
density regions where star formation is currently occur- 
ring were formed through the interaction of newly formed 
stars and the surrounding molecular material. Carpenter 
et al. (1995a) split the cloud into 7 separate regions and 
find that 3 regions (S247, S252, and S254-258) contain 
active star formation based on higher H2 column densi- 
ties (> 10^^ cm~^) and warmer kinetic temperatures (10 
to 30 K). 

Carpenter et al. (1995b) identified 11 clumps within 
the ^^CO filaments using CS(J = 2 — 1) observations 
(55" resolution), all with masses greater than 100 Mq. 
They found that at least eight of the 1 1 clumps have asso- 
ciated IRAS point sources. Carpenter et al. (1995b) also 
conducted a near-infrared imaging survey of a 30' x 45' 
area surrounding S247, including three of the CS clumps. 
They found that each of the three clumps contains a clus- 
ter of stars in the near-infrared images. 

Chavarria et al. (2008) identified 510 young stellar ob- 
jects (YSOs) with near or mid-IR excess in the S254-258 



3.1.1. Observations 
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Fig. 1. — The Gemini OBI molecular cloud. Extinction maps of 
Dobashi et al. (2005) are shown in inverted grayscale ranging from 
-0.5 to 4.0 magnitudes. The OB association members (Humphreys 
1978) are shown as stars. The HII regions (Sharpless 1959) and IC 
443 are shown as circles sized to match the extent of the region. 
The dashed line marks the area mapped in ^^CO and ^^CO by 
Carpenter et al. (1995a), and the solid line marks the extent of the 
1.1 mm BGPS maps. 



complex using a combination of Spitzer IRAC and near- 
infrared data. They classify 87 Class I sources, and 165 
Class II sources, and find that 80% of the YSOs are in 
clusters surrounded by a more evenly distributed popu- 
lation. 

Figure [1] shows the extinction maps of Dobashi et al. 
(2005) for the Gemini OBI molecular cloud. OB associa- 
tion members are shown as stars (Humphreys 1978), HII 
regions (Sharpless 1959) and IC 443 are shown as circles 
with sizes corresponding to the source size. The extent 
of the i^CO and "CO maps of Carpenter et al. (1995a) 
is marked by the dashed line, while the solid line shows 
the extent of the BGPS map presented in this paper. 

3. SURVEYS 

3.1. Bolocam 1.1 mm Galactic Plane Survey 

The BGPS^"^ has surveyed approximately 170 square 
degrees of the Galactic Plane using the Bolocam in- 
strument on the Caltech Submillimeter Observatory^'* 
(CSO) on Mauna Kea. Aguirre et al. (2010) present 
specific details regarding the BGPS and methods and 
Rosolowsky et al. (2010) detail the BGPS catalog. The 
images and catalog have been made public and are hosted 
by the Infrared Processing and Analysis Center via the 
NASA/IPAC Infrared Science Archive" (IPAC). 



^^ See http://milkyway.colorado.edu/bgps/ 

'^^ The Caltech Submillimeter Observatory is supported by the 



NSF. 
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See |http://irsa.ipac.caltech.edu/data/BOLOCAM_GPS/| 



The Gem OBI region was observed in two 3°xl° 
groups of observations, centered at {l,h) = (189°, 0.5°) 
and (192°, 0°), with a shallow observation centered be- 
tween the two main groups to aid in mosaicing the im- 
ages (see Figure [21) . The centers were selected to place 
the most active star-forming regions (S247, S252, and 
S254-258) within our mapped region. Observations were 
carried out at the Caltech Submillimeter Observatory 
(CSO) as part of the BGPS between 9 and 13 September 
2007. The survey used the Bolocam instrument*^, which 
has a hexagonal array of 144 bolometers. The Bolocam 
filter is centered at 268 GHz with a bandwidth of 45 GHz, 
and was designed to exclude the *^C0(2 — >■ 1) emission 
line. Each bolometer beam is well approximated by a 
Gaussian with a FWHM of 31", and the bolometer sep- 
aration results in an instantaneous field of view of 7.5'. 
For the BGPS, the effective beam size after processing 
many scans is 33" for the current release (Aguirre et al. 
2010). For more instrument specific details, see Glenn et 
al. (2003) and Haig et al. (2004). Observations were ob- 
tained using the raster scan mode, with alternating scans 
directed along I and h. The scan speed was 120" s~*, and 
chopping was not utilized in these observations, thereby 
retaining some sensitivity to large scale structure up to 
the angular size of the array (see Aguirre et al. 2010 for 
more detailed analysis). 

3.1.2. Reduction 

At millimeter wavelengths, the detected signal is dom- 
inated by the atmosphere, making the sky cleaning 
algorithm extremely important. Bolocam's individual 
bolometer beams overlap substantially in the lower at- 
mosphere such that they probe similar columns of at- 
mosphere, facilitating sky subtraction. The iterative 
mapping method utilizes principal component analysis 
(PCA) and considers all observations for a given field si- 
multaneously. In our iterative mapping method, the tem- 
porally correlated signal is first subtracted from the time 
streams. This subtraction removes both the atmosphere 
and all astrophysical structure larger than the array field 
of view. The remaining time stream, which in principle 
should contain only astrophysical signal, is then mapped 
with 7.2" X 7.2" pixels, deconvolved from the beam, and 
returned to a time stream. The time stream containing 
only astrophysical signal is subtracted from the original 
to produce a third time stream (summing the subtracted 
astrophysical signals to make a source map) that repre- 
sents an estimate of the noise. The above procedure is 
then repeated for 50 iterations on the deconvolved as- 
trophysical signal time stream to produce the final map. 
This procedure can be fine-tuned to balance the retention 
of large-scale structure against production of artifacts 
(such as negative bowls) around bright sources. The it- 
erative mapping method resulted in a 1-a RMS noise of 
0.070 Jy/beam in the Gemini OBI image. For more spe- 
cific information on the BGPS reduction, see Aguirre et 
al. (2010). 

3.1.3. Images 
^^ See |http://www. cso.caltech.edu/bolocam /I 
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Fig. 2. — The BGPS 1.1 mm continuum image of the Gemini OBI molecular cloud shown in an inverted, linear grayscale ranging 
-0.2 to 0.75 Jy beani"^. Groups 1, 2, and 3 are shown in more detail in following figures. 

The iteratively mapped BGPS image is shown in Fig- 
ure[2] The milhmeter sources in our maps are highly clus- 
tered and are found in three groups, while the remainder 
of the image is strikingly devoid of millimeter emission. 
This clustering is in contrast to the inner Galaxy regions 
also included in the BGPS, where a similarly sized im- 
age would contain hundreds of sources [see Aguirre et 
al. (2010) and J. Bally et al. (2010)]. The three groups 
correspond to S247 (referred to as Group 1; Figure |3]), 
S252 (referred to as Group 2; Figure S]), and S254-258 
(referred to as Group 3; Figure [5]). 

3.2. NHg Survey of BGPS Sources 

The NH3 pointings are displayed in Figures [3l-l5] as 
crosses, diamonds, and asterisks. Since the NH3 point- 
ings were chosen prior to completion of the BGPS cata- 
log, they do not necessarily correspond to the peak posi- 
tions of the millimeter emission determined by the source 
extraction algorithm (see §4.1|) . Crosses mark the posi- 
tions closest to the peak of the millimeter emission, dia- 
monds mark extra NH3 pointings within each source, and 
the asterisks mark the positions of NH3 pointings which 
correspond to a mm source but do not have a detection 
in any of the NH3 transitions. 

3.2.1. Observations 

We observed 56 NH3 pointings in the Gem OBI re- 
gion using the Robert F. Byrd Green Bank Telescopc"'^^ 

^^ The GET is operated by the National Radio Astronomy Ob- 
servatory, which is a facility of the National Science Foundation, 
operated under cooperative agreement by Associated Universities, 
Inc. 
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Fig. 3. — The BGPS image of the region marked Group 1 in Fig- 
ure [2] in inverted grayscale ranging from -0.2 to 0.75 Jy beam~ . 
The BGPS sources with NH3 spectra are shown as black ellipses, 
while BGPS sources without NH3 spectra are shown as white el- 
lipses. The NH3 pointings are shown as crosses and diamonds. 
Crosses denote pointings closest to the peak of the millimeter emis- 
sion, while diamonds denote additional NH3 pointings within each 
BGPS source. 
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Fig. 4. — The BGPS image of the region marked Group 2 in Fig- 
ure [2] in inverted grayscale ranging from -0.2 to 0.75 Jy beam . 
The BGPS sources with NH3 spectra are shown as black ellipses, 
while BGPS sources without NH3 spectra are shown as white el- 
lipses. The NH3 pointings are shown as crosses, diamonds and 
asterisks. Crosses denote pointings closest to the peak of the mil- 
limeter emission, diamonds denote additional NH3 pointings within 
each BGPS source, and asterisks denote non-detections. Source el- 
lipses can overlap since these are the moment approximations to 
source sizes rather than the boundaries established by the extrac- 
tion algorithm. 



(GBT) on 14 and 15 February 2008 in two separate 
blocks totaling 3.6 hours. The 23 GHz zenith opacity 
for the Feb. 14 observations was T23 ^ 0.04; and for 
Feb. 15, T23 ^ 0.11. The observations were broken 
into two minute integrations of individual targets. We 
used Feed 1 of the high-frequency K-band receiver on the 
GBT as the front end for the observations. The GBT 
spectrometer served as the back end, which we config- 
ured to observe both polarizations of the feed with 50 
MHz bandwidth for four separate IFs. The total band- 
width was sampled in 8192 channels for a 6.1 kHz (^0.08 
km s""'^) channel width. We used two different sets of 
IFs. All scans contained the NH3(1,1) (23.6944955(1) 
GHz) and NH3(2,2) (23.7226333(1) GHz) lines. Obser- 
vations on February 14 (10 pointings) also included the 
C2S (Nj = 2i ^ lo, 22.344033(1) GHz, Yamamoto et al. 
1990) and NH3(4,4) (24.1394169(1) GHz) lines. All NH3 
frequencies were taken from Lovas & Dragoset (2003). 
There were no significant detections of the C2S line and 
only one marginal detection of the NH3(4,4) line, so the 
spectrometer configuration was changed for the follow- 
ing days' observations. The observations on February 
15 included the H2O (22.344033 GHz), and NH3(3,3) 
(23.8701279(5) GHz) transitions in the remaining two 
IFs. Ten pointings were observed on 14 February 2008 
with C2S and NH3(4,4), and all 56 pointings were ob- 
served with H2O and NH3(1,1), (2,2) and (3,3). Ten 
pointings were observed in all six lines (four NH3 lines. 
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Fig. 5. — The BGPS image of the region marked Group 3 in Fig- 
ure [2] in inverted grayscale ranging from -0.2 to 0.75 Jy beam . 
The BGPS sources with NH3 spectra are shown as black ellipses, 
while BGPS sources without NH3 spectra are shown as white el- 
lipses. The NH3 pointings are shown as crosses, diamonds and 
asterisks. Crosses denote pointings closest to the peak of the mil- 
limeter emission, diamonds denote additional NH3 pointings within 
each BGPS source, and asterisks denote non-detections. Source el- 
lipses can overlap since these are the moment approximations to 
source sizes rather than the boundaries established by the extrac- 
tion algorithm. 

H2O, and C2S). We observed each source with a sym- 
metric, in-band frequency switch (5 MHz throw). 

Observations of the Gem OBI region were conducted 
immediately following observations for a pointing model 
correction on each night. We calibrated the data with 
the injection of a noise signal periodically throughout 
the observations. Because of slow variations in the power 
output of the noise diodes and their coupling to the sig- 
nal path, we measured the strength of the noise signal 
through observations of a source with known flux (the 
NRAO flux calibrator 3C 48, S^ = 1.13 Jy). We re- 
peated the flux calibration observations during every ob- 
serving session to detect any changes in the calibration 
sources, finding no significant variations over the course 
of our run (< 10%). The beam size of the GBT is 31" at 
these frequencies, which projects to a si ze of 0.32 pc at 
the distance of Gem OBI (2.10 kpc, see %.![ Reid et al. 
(2009)). The GBT NH3 observations are weU matched 
to the BGPS observations since the GBT beam size is 
only a few arcseconds smaller than the effective BGPS 
beam size. 

3.2.2. Data Reduction 

The ammonia observations were reduced using the 
GBTIDL reduction package. We extracted each scan, 
using the package to fold the spectrum to remove the ef- 
fects of the frequency switching and establish the data 
on the Ta scale by calibrating the noise diodes. We sub- 
tracted a linear baseline from the spectrum from each IF 



TABLE 1 

Mean RMS for Spectral Line 

Observations 





mean RMS 


stdev RMS 


Line 


(mK) 


(mK) 


NH3{1,1) 


82 


21 


NH3(2,2) 


81 


21 


NH3(3,3) 


97 


11 


NH3(4,4) 


52 


1 


H2O 


172 


26 


C2S 


67 


4 



with the exception of the IF containing the water and 
C2S Hnes where a second order basehne was required. 
We then scaled to the T^ scale using estimates of the 
atmospheric opacity interpolated in the frequency and 
time domains at 22 and 23 GHz from models of the at- 
mosphere derived using weather data.^^ To reach the 
Tmb scale, the spectra are divided by the main beam 
efficiency of the GET, which is rjmb ~ 0.81 at these fre- 
quencies. We averaged together both polarizations; if a 
source was observed multiple times, each observation was 
calibrated separately and the resulting spectra were then 
averaged together. The mean and standard deviation of 
the RMS noise per channel for each of the 6 observed 
lines are listed in Table [TJ 

4. PARAMETER ESTIMATION 

4.1. BGPS Parameter Estimation 

The BGPS millimeter source properties are extracted 
from the iterative maps using the BGPS source extrac- 
tion software, Bolocat (Rosolowsky et al. 2009). Bolo- 
cat utilizes a seeded watershed method (also referred to 
as marker-controlled segmentation; Soille 1999) which is 
comprised of two main steps. First, regions with emis- 
sion above a signal-to-noise ratio of 2.5 are identified 
and expanded through a nearest-neighbor algorithm to 
include adjacent, lower signal-to-noise pixels. Second, 
each identified region is subdivided into further sub- 
structures based on the contrast between local maxima. 
Bolocat then determines source properties from emission- 
weighted moments. Two sets of coordinates are deter- 
mined for each source: peak and centroid coordinates. 
The peak coordinates correspond to the peak of the 1.1 
mm emission and are determined from median smoothed 
maps. The centroid coordinates are given by the first 
emission-weighted moments. The major and minor axis 
dispersions {(Jmaj and a mm) are given by the second mo- 
ments. The object radius (Robj) is calculated from 
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where abeam = Qfwhm/VS In 2 and r? is a factor relat- 
ing the axis dispersions to the true size of the source, 
found to be 2.4. The value ry — 2.4 is adopted as the 
median value derived from measuring the observed axis 
dispersions compared with the true radius for a variety of 
simulated emission profiles spanning a range of density 
distributions, sizes relative to the beam and signal-to- 
noise ratios. The ellipses in Figures [3HS] mark the po- 
sitions and sizes of the extracted BGPS sources. White 
ellipses denote BGPS sources which do not contain any 

^® See |http : / /www. gb .nrao.edu/ ~rmaddale/ Weather | 



TABLE 2 

Aperture Flux 
Corrections 



Aperture 
(") 



S^(int)/S'^(ap) 



40 
80 

120 



1.46 
1.04 
1.01 



NH3 pointings or only contain a pointing which was not 
detected, and black ellipses denote BGPS sources which 
contain at least one detected NH3 pointing. 

Source fiux densities are determined by aperture pho- 
tometry within standard apertures, 40", 80", 120", and 
2Robj diameter apertures denoted as 5*1,(40"), 5,^(80"), 
5*1/(120"), and Si,{ohi). An integrated fiux density, 
5i/(int), is also reported for each source and is the sum 
of all pixels assigned to a given source by the source ex- 
traction algorithm divided by the number of pixels per 
beam. The intensity in the BGPS maps was calibrated 
based on the observed properties of planets. Aguirre et 
al. (2010) have compared the BGPS data to other mm 
continuum surveys and found that we must multiply our 
flux densities by 1.5±0.15 in order for our fiux densities 
to agree with the other surveys (see Aguirre et al. (2010) 
for more details). In this paper we apply the flux calibra- 
tion correction of 1.5±0.15 to the vl.O BGPS products 
presented in Rosolowsky et al. (2010) and released on 
the IPAC website. All data products and properties pre- 
sented in this paper include this flux calibration factor. 

Aperture corrections must be applied to the photome- 
try in order to account for power at large radii due to the 
non-gaussian side lobes of the CSO beam. The aperture 
corrections were determined by comparing the integrated 
flux density for known point sources to the calculated 
aperture fiux densities. The aperture correction values 
are listed in Table [H As can be seen from the aperture 
corrections, the 40" aperture underestimates the flux 
density by almost 50%. The 80" aperture only under- 
estimates the flux density by approximately 4%, and the 
120" aperture requires only a 1% correction. The aper- 
tures corresponding to 2Robj require varying aperture 
corrections, while 5,y(int) represents all fiux in the 1.1 
mm maps down to the la level and does not require an 
aperture correction factor. The aperture corrections were 
determined based on sources very small compared to the 
BGPS beam and therefore will slightly overestimate the 
flux densities of sources which are extended compared 
to the beam size. For more details regarding source ex- 
traction and parameter estimation, see Rosolowsky et al. 
(2010). 

4.2. NH3 Parameter Estimation 

We estimate physical properties of the molecular gas by 
fitting a simple model to each calibrated spectrum, fol- 
lowing the method described in Rosolowsky et al. (2008; 
hereafter R08). The model estimates the ammonia emis- 
sion that would be observed from a slab of molecular gas 
with a fixed optical depth in the NIl3(l,l) line of r, LSR 
velocity i'lsR; velocity dispersion ay, radiative excitation 
temperature Tg^-, and kinetic temperature Tr-. For clar- 
ity, we note that r here represents the total opacity in 
the (1,1) transition which is then apportioned into the 18 
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Fig. 6. — The distribution of extracted source radii for all 55 
BGPS sources. The white histogram marks all 55 sources while 
the gray histogram marks the unresolved sources for which we have 
set the radius to an upper-limit of the beam size. The dotted lines 
mark the radii of the aperture sizes used in the BGPS catalog, 
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hyperfine components of the line according to their sta- 
tistical weights. The central cluster of ammonia hyper- 
fine components, which is commonly studied in isolation, 
contains half the total optical depth, T„iain — 0.5t. 

The model spectrum provides a simultaneous estimate 
of all the ammonia parameters; however, specific aspects 
of the spectral data control each derived quantity. For 
example, the ratio of strengths of the (1,1) and (2,2) lines 
yields an estimate of the kinetic temperature. The model 
accounts for the line widths by assuming a Gaussian ve- 
locity distribution in the underlying slab of material so 
that the model spectrum accounts for line broadening by 
hyperfine structure and optical depth efi^ects. The differ- 
ence between the line strength on the Tmb scale and the 
kinetic temperature provides information on the radia- 
tive excitation of the line Tex under the assumption that 
the emission fills the telescope beam (see ^ 36. 61 for more 
details). As noted in R08, the parameters T^x and r are 
strongly covariant and cannot be determined separately 
in the low optical depth limit. In these cases, only the 
product tT^x is determined. At higher optical depths, 
the relative strengths of the central and satellite hyper- 
fine transitions in the (1,1) line provide an independent 
estimate of the optical depth. 

Using these parameters as a description of the under- 
lying molecular gas, a model spectrum can be generated 
that adequately describes emission we observe from the 
Gem OBI sources including the hyperfine structure of the 
ammonia transitions. The best-fit set of parameters is 
determined by a least-squares optimization of the model 
spectrum to fit to all observed ammonia transitions si- 
multaneously. The primary advantage of this approach is 
the straight-forward estimation of errors and covariances 
in the fit parameters. 

The spectral model has two significant differences from 
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Fig. 7. — The distribution of distances between the 55 BGPS 
sources out to 500". The dotted lines mark the diameters of the 
apertures used in the BGPS catalog, Bap = 40", 80", and 120". 
Only approximately 5 sources will suffer from contamination with 
the 120" aperture. 

that of the R08 study. In particular, the inclusion of 
NIl3(3,3) lines in the analysis necessitates treating the 
ortho- and para-NHa states in more detail. Rather than 
simply adopting a full thermal equilibrium as done by 
R08, we treat the two sets of states as different species 
with an ortho-to-para ratio of 1:1, which is the high tem- 
perature limit [Tk > 40 K; Takano et al. 2002). Given 
that the kinetic temperature derived from only the (1,1) 
and (2,2) lines provides a good estimate of the (3,3) line 
in most cases, this assumption is adequate. As such, 
partition functions for the two species are decoupled and 
contain only the ortho- and the para- states: 



Zn = 



l+^2(2J+l)exp|- 

7 ^ 



-h[BJ{J + l) + {C-B)J^] 



for J = 3,6,9, 

Zp = 
^(2J-}-l)exp 



fcTfe 



-h[BJ{J+l) + {C-B)J^] 



(2) 



kTk 



for J = 1,2,4,5,, 



(3) 



where we have only included the metastable states given 
hy J = K. The lifetimes of the non-metastable states 
are relatively short and the levels are not well populated. 
The partition functions are subject to the constraint that 
Zo — Zp^ representing our assumption of the 1:1 ortho- 
to-para ratios. In addition to a different partition func- 
tion, we also model the effects of the frequency switch- 
ing explicitly in the model spectrum since the frequency 
switch is smaller than in the ROB study. We build a fre- 
quency switched version of the spectrum by shifting it 



TABLE 3 

Millimeter Sources 



ID 




BGPS ID 


Group 




Number 


Source 


Number 


Number 


Comments 


01 


G189.646+00.131 


7469 


2 




02 


G189.659+00.185 


7470 


2 




03 


G189.682+00.185 


7471 


2 


IRAS 060474-2040 


04 


G189.782+00.265 


7475 


2 


IRAS 060534-2036 


05 


G189.116+00.643 


7467 


1 


IRAS 06052-1-2122 


06 


G189.788+00.281 


7479 


2 




07 


G189.713+00.335 


7472 


2 




08 


G189.789+00.291 


7480 


2 




09 


G189.744+00.335 


7473 


2 




10 


G189.782+00.323 


7476 


2 




11 


G189.836+00.303 


7485 


2 


IRAS 060554-2034 


12 


G189.950+00.231 


7491 


2 




13 


G189.776+00.343 


7474 


2 


S252A, IRAS 06055-1-2039 


14 


G189.834+00.317 


7484 


2 




15 


G189.888+00.303 


7489 


2 


IRAS 06056-f 2032 


16 


G189.804+00.355 


7481 


2 




17 


G189.030+00.781 


7465 


1 


IRAS 060564-2131 


18 


G189.831+00.343 


7843 


2 


J06084309+2036182 


19 


G189.879+00.319 


7487 


2 




20 


G189.885+00.319 


7488 


2 




21 


G189.810+00.369 


7482 


2 




22 


G189.032+00.793 


7466 


1 


IRAS 060564-2131 


23 


G189.921+00.331 


7490 


2 




24 


G189.015+00.823 


7464 


1 




25 


G188.991+00.859 


7463 


1 




26 


G188.948+00.883 


7461 


1 


S252 


27 


G189.951+00.331 


7492 


2 




28 


G189.783+00.433 


7477 


2 


IRAS 060594-2042 


29 


G188.975+00.911 


7462 


1 




30 


G189.990+00.353 


7493 


2 




31 


G189.783+00.465 


7478 


2 




32 


G188.792+01.027 


7460 


1 


IRAS 060614-2151 


33 


G190.006+00.361 


7494 


2 


IRAS 060614-2028 


34 


G189.864+00.499 


7486 


2 


IRAS 060634-2040 


35 


G189.231+00.893 


7468 


1 


IRAS 060654-2124 


36 


G190.054+00.533 


7496 


2 


IRAS 060684-2030 


37 


G190.044+00.543 


7495 


2 




38 


G190.063+00.679 


7497 


2 




39 


G190.192+00.719 


7499 


2 




40 


G190.171+00.733 


7498 


2 


IRAS 060784-2030 


41 


G190.240+00.911 


7500 


2 


IRAS 060874-2031 


42 


G192. 629-00.157 


7505 


3 




43 


G192. 602-00.143 


7503 


3 


S256, IRAS 060964-1757 


44 


G192. 629-00.123 


7504 


3 




45 


G192. 581-00.043 


7501 


3 


S255N 


46 


G192. 662-00.083 


7507 


3 




47 


G192. 596-00.051 


7502 


3 


S255IR, J061253304-1759215 


48 


G192.644+00.003 


7506 


3 


S258 


49 


G192.719+00.043 


7508 


3 


IRAS 061054-1756 


50 


G192.764+00.101 


7509 


3 




51 


G192.816+00.127 


7510 


3 




52 


G192.968+00.093 


7511 


3 




53 


G193.006+00.115 


7514 


3 




54 


G192.981+00.149 


7512 


3 


J061423594-1744469 


55 


G192.985+00.177 


7513 


3 





±5 MHz in frequency and subtracting half of the aniph- 
tude in the model. Thus, we fit both the peaks and the 
negative features in the observed spectra, resulting in a 
fit which is more sensitive to the low-intensity emission. 
Emission from the H2O and C2S lines is not treated in 
detail by the model but rather by simply calculating its 
velocity-integrated intensity. Since the model is neces- 
sarily simple, complicated sources will be represented by 
average properties. However, the model produces good 
fits for the wide variety of sources found in the observa- 
tions. 



5. RESULTS 



5.1. 1.1 mm Results 

We have extracted 55 millimeter continuum sources 
from the BGPS maps. The 1.1 mm sources are listed 
in Table |31 along with their ID number in the complete 
BGPS source catalog, their corresponding group number 
(see §3.1.3p . and any coincident IRAS sources, 2MASS 
sources and well known millimeter sources. Table 2] lists 
extracted 1.1 mm source properties including a running 
source ID number (column 1), the source name (based on 
the Galactic coordinates of the peak emission; Column 
2), centroid RA and DEC (columns 3 and 4), RA and 
DEC of the peak emission (columns 5 and 6), size of semi- 
major and semi-minor axes in arcseconds (columns 7 and 



TABLE 4 
Observed 1.1 mm Properties 



ID 




RA 


Dec 


peak RA 


peak Dec 


^major 


-tv'Tninor 


PA 


Rob/ 


S,.(120") 


S^(int) 


number 


Source 


{J2000) 


(J2000) 


(J2000) 


{J2000) 


(") 


(") 


(°) 


(pc) 


(mjy) 


(mJy) 


01 


G189.646+00.131 


06 07 30.6 


20 40 00.6 


06 07 30.9 


20 39 46.1 


33.0 


19.9 


74 


0.50 


936 


(250) 


988 (120) 


02 


G189.659+00.185 


06 07 44.4 


20 40 49.6 


06 07 44.7 


20 40 36.6 


20.5 


14.1 


139 


< 0.17 


1010 


(270) 


370 (91) 


03 


G189. 682+00. 185 


06 07 46.9 


20 39 39.6 


06 07 47.5 


20 39 27.4 


33.5 


19.2 


138 


0.49 


1509 


(300) 


1700 (150) 


04 


G189. 782+00. 265 


06 08 17.3 


20 36 37.9 


06 08 17.9 


20 36 32.6 


28.2 


18.0 


145 


0.40 


1380 


(300) 


1190 (140) 


05 


G189. 116+00. 643 


06 08 19.9 


21 22 26.3 


06 08 19.8 


21 22 29.2 


31.2 


19.9 


4 


0.49 


1190 


(320) 


1090 (130) 


06 


G189. 788+00. 281 


06 08 22.8 


20 36 38.0 


06 08 22.2 


20 36 41.7 


21.7 


17.6 


43 


0.32 


1630 


(330) 


1019 (140) 


07 


G189. 713+00. 335 


06 08 25.4 


20 42 04.8 


06 08 25.1 


20 42 08.9 


26.3 


20.6 


158 


0.45 


1360 


(300) 


1160 (140) 


08 


G189. 789+00. 291 


06 08 25.8 


20 37 04.9 


06 08 24.7 


20 36 52.8 


16.5 


9.1 


13 


< 0.17 


1430 


(330) 


347 (94) 


09 


G189. 744+00. 335 


06 08 29.0 


20 41 04.7 


06 08 28.8 


20 40 34.5 


28.6 


20.9 


72 


0.48 


2960 


(420) 


1960 (210) 


10 


G189. 782+00. 323 


06 08 29.3 


20 38 13.1 


06 08 30.9 


20 38 13.9 


28.5 


22.7 


38 


0.51 


3190 


(450) 


2240 (220) 


11 


G189. 836+00. 303 


06 08 30.5 


20 34 33.3 


06 08 33.1 


20 34 49.0 


34.2 


23.4 


142 


0.59 


2480 


(380) 


1780 (210) 


12 


G189. 950+00. 231 


06 08 31.0 


20 26 49.9 


06 08 31.1 


20 26 44.4 


34.4 


22.6 


17 


0.58 


1670 


(310) 


1910 (160) 


13 


G189. 776+00. 343 


06 08 34.4 


20 39 18.2 


06 08 34.6 


20 39 07.7 


38.5 


26.8 


151 


0.70 


8930 


(850) 


11100 (410) 


14 


G189. 834+00. 317 


06 08 35.0 


20 35 10.5 


06 08 36.0 


20 35 19.7 


23.4 


17.1 


32 


0.33 


2730 


(400) 


1530 (170) 


15 


G189. 888+00. 303 


06 08 38.9 


20 32 06.5 


06 08 39.6 


20 32 05.3 


20.6 


13.5 


165 


<0.17 


1370 


(290) 


537 (110) 


16 


G189. 804+00. 355 


06 08 39.7 


20 38 02.0 


06 08 40.8 


20 38 00.5 


23.4 


22.2 


40 


0.44 


8060 


(790) 


4700 (320) 


17 


G189. 030+00. 781 


06 08 40.2 


21 31 00.0 


06 08 40.1 


21 31 00.8 


26.9 


22.8 


85 


0.50 


11100 


(990) 


11200 (410) 


18 


G189. 831+00. 343 


06 08 41.3 


20 36 04.8 


06 08 41.6 


20 36 11.4 


36.2 


28.9 


73 


0.71 


5780 


(610) 


6420 (350) 


19 


G189. 879+00. 319 


06 08 43.2 


20 33 49.9 


06 08 42.2 


20 32 58.4 


39.2 


21.9 


58 


0.61 


2390 


(350) 


2400 (210) 


20 


G189. 885+00. 319 


06 08 44.5 


20 32 12.0 


06 08 42.9 


20 32 39.5 


32.5 


22.7 


1 


0.56 


2450 


(360) 


2080 (200) 


21 


G189. 810+00. 369 


06 08 44.6 


20 38 18.5 


06 08 44.7 


20 38 06.1 


54.4 


31.8 


66 


0.95 


8090 


(780) 


11300 (480) 


22 


G189. 032+00. 793 


06 08 45.7 


21 31 51.7 


06 08 43.1 


21 31 15.4 


30.9 


23.7 


39 


0.56 


11100 (1000) 


6690 (430) 


23 


G189. 921+00. 331 


06 08 49.2 


20 31 46.8 


06 08 50.1 


20 31 07.1 


25.0 


14.9 


21 


0.25 


2500 


(380) 


697 (160) 


24 


G189. 015+00. 823 


06 08 49.7 


21 33 49.2 


06 08 47.9 


21 32 58.1 


34.9 


16.2 


67 


0.39 


2690 


(390) 


1660 (190) 


25 


G188. 991+00. 859 


06 08 53.1 


21 35 32.0 


06 08 53.0 


21 35 16.5 


30.2 


12.6 


51 


< 0.17 


977 


(270) 


707 (110) 


26 


G188. 948+00. 883 


06 08 53.4 


21 38 23.8 


06 08 52.9 


21 38 16.9 


29.0 


24.4 


89 


0.55 


12899 ( 


1200) 


14599 (450) 


27 


G189. 951+00. 331 


06 08 53.6 


20 29 37.3 


06 08 53.8 


20 29 32.6 


53.5 


34.9 


157 


0.99 


5980 


(620) 


10800 (390) 


28 


G189. 783+00. 433 


06 08 55.5 


20 41 16.0 


06 08 55.8 


20 41 19.6 


43.6 


22.4 


44 


0.66 


1439 


(300) 


1840 (160) 


29 


G188. 975+00. 911 


06 09 03.6 


21 38 03.8 


06 09 02.7 


21 37 37.5 


26.6 


16.2 


56 


0.33 


718 


(260) 


703 (95) 


30 


G189. 990+00. 353 


06 09 03.7 


20 28 18.0 


06 09 03.4 


20 28 11.3 


22.3 


14.0 


3 


<0.17 


886 


(280) 


499 (91) 


31 


G189. 783+00. 465 


06 09 04.4 


20 42 12.6 


06 09 03.0 


20 42 15.4 


30.0 


19.8 


124 


0.47 


917 


(270) 


949 (120) 


32 


G188. 792+01. 027 


06 09 06.7 


21 50 48.5 


06 09 06.0 


21 50 39.1 


27.2 


22.2 


83 


0.49 


4240 


(570) 


4650 (250) 


33 


G190. 006+00. 361 


06 09 08.4 


20 27 33.7 


06 09 07.2 


20 27 34.9 


30.5 


20.1 


110 


0.48 


887 


(260) 


760 (110) 


34 


G189. 864+00. 499 


06 09 20.7 


20 39 30.8 


06 09 20.6 


20 39 02.7 


50.2 


36.5 


27 


0.99 


3850 


(450) 


6969 (300) 


35 


G189. 231+00. 893 


06 09 30.4 


21 23 44.7 


06 09 30.6 


21 23 39.8 


17.4 


15.4 


100 


0.20 


464 


(260) 


427 (79) 


36 


G190. 054+00. 533 


06 09 51.1 


20 30 14.5 


06 09 51.7 


20 30 03.2 


52.7 


23.4 


21 


0.75 


2220 


(340) 


2930 (220) 


37 


G190. 044+00. 543 


06 09 54.6 


20 30 56.2 


06 09 52.7 


20 30 52.2 


25.1 


17.7 


22 


0.37 


2420 


(360) 


1019 (170) 


38 


G190. 063+00. 679 


06 10 27.0 


20 33 48.5 


06 10 25.7 


20 33 45.6 


20.3 


12.9 


149 


<0.17 


418 


(230) 


229 (64) 


39 


G190. 192+00. 719 


06 10 50.3 


20 27 51.6 


06 10 50.5 


20 28 11.4 


16.4 


12.1 


29 


<0.17 


685 


(260) 


316 (74) 


40 


G190. 171+00. 733 


06 10 51.3 


20 29 57.7 


06 10 51.2 


20 29 38.8 


45.1 


16.9 


32 


0.49 


1540 


(300) 


2090 (160) 


41 


G190. 240+00. 911 


06 11 40.4 


20 31 26.7 


06 11 39.6 


20 31 12.8 


28.0 


17.2 


73 


0.38 


726 


(260) 


626 (98) 


42 


G192. 629-00.157 


06 12 33.6 


17 54 48.7 


06 12 33.5 


17 54 40.4 


21.9 


16.7 


57 


0.30 


796 


(290) 


558 (96) 


43 


G192. 602-00.143 


06 12 33.6 


17 56 19.8 


06 12 33.2 


17 56 33.1 


25.9 


16.8 


11 


0.35 


894 


(300) 


868 (110) 


44 


G192. 629-00.123 


06 12 39.8 


17 55 28.8 


06 12 41.0 


17 55 39.2 


35.6 


21.3 


16 


0.56 


709 


(290) 


740 (110) 


45 


G192. 581-00.043 


06 12 52.8 


18 00 32.6 


06 12 52.9 


18 00 29.0 


25.2 


21.9 


174 


0.46 


14600 ( 


1300) 


11400 (450) 


46 


G192. 662-00.083 


06 12 53.1 


17 55 21.0 


06 12 53.7 


17 55 07.3 


33.5 


11.7 


43 


<0.17 


565 


(310) 


463 (85) 


47 


G192. 596-00.051 


06 12 54.1 


17 58 54.1 


06 12 52.8 


17 59 31.0 


46.5 


29.5 


22 


0.83 


15400 (1400) 


17300 (590) 


48 


G192. 644+00. 003 


06 13 09.8 


17 58 40.3 


06 13 10.6 


17 58 32.6 


28.1 


22.8 


45 


0.51 


1330 


(320) 


1380 (140) 


49 


G192. 719+00. 043 


06 13 28.0 


17 56 02.1 


06 13 28.6 


17 55 41.6 


43.2 


27.3 


120 


0.76 


2320 


(390) 


3580 (220) 


50 


G192. 764+00. 101 


06 13 46.7 


17 55 00.4 


06 13 46.8 


17 55 02.5 


39.6 


21.0 


146 


0.59 


2020 


(370) 


2480 (180) 


51 


G192. 816+00. 127 


06 13 58.8 


17 52 51.2 


06 13 58.8 


17 53 03.0 


20.1 


13.2 


39 


<0.17 


615 


(320) 


505 (84) 


52 


G192. 968+00. 093 


06 14 08.3 


17 44 04.6 


06 14 09.6 


17 44 03.9 


25.4 


16.2 


166 


0.32 


682 


(310) 


620 (93) 


53 


G193. 006+00. 115 


06 14 18.8 


17 42 36.9 


06 14 19.1 
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17 44 56.1 


49.6 


25.5 
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06 14 30.2 


17 45 34.0 


06 14 30.4 


17 45 31.6 


18.4 


14.3 


85 


< 0.17 


844 


(320) 


516 (90) 



Note. — Errors arc given in parcntlicscs. 

Sources tliat arc unresolved with tlie BGPS beam are assigned the beam size as an upper-limit on the object radius. 



8), position angle of best- fit ellipse measured counter- 
clockwise from the x-axis {i = 0°) in Galactic coordi- 
nates (column 9), the deconvolved source radius Robj 
(column 10), and 5*^(120") and S^{mt) (columns 11 and 
12). Some sources are sufficiently small that a decon- 
volved radius cannot be determined, and we assume a 
source diameter equal to the beam FWHM as an upper 
limit on the radius for these objects. 

5'y(int) and Robj adequately characterize the sources 
but these two quantities are highly dependent upon the 



source extraction algorithm and are not easily compared 
with other observations. In order to support comparison 
of our results with other studies, we will also present 
a flux within a well-characterized and commonly used 
aperture. 

To evaluate which aperture most closely represents the 
true mm flux, we consider the derived object radii as 
well as the distance between all BGPS sources. Fig- 
ure ini shows the distribution of derived source radii for 
the 55 BGPS sources (white histogram) and the unre- 
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Fig. 8.— H20,NH3(1,1), (2,2), and (3,3) spectra for G189. 776-1-0.343 (ID #13). Top panel: Black is observed NH3(1,1) spectrum and 
green is the model fit. From left to right, the bottom panels show H2O, NH3(2,2) and NH3(3,3) spectra. Figures 8b-at are available in the 
online version of the Journal. 



solved sources where we have assigned the beam size as 
an upper- limit for the source radius (gray histogram). 
The dotted hues mark the radii of the apertures in which 
we present photometry in the BGPS catalog. There is 
a large spread in the object radii, with the largest being 
almost 100". Figure [7] shows the distribution of angular 
separation between the emission maxima of all 55 BGPS 
sources out to 500". The dotted lines mark the diameters 
of our chosen apertures. 

To characterize the mm flux well, we should choose 
an aperture size greater than the largest source radius 
yet small enough to reduce the number of sources which 
will include flux from their neighbors. If we choose an 
aperture much larger than the source and no other mm 
sources fall within that aperture, then the flux will still be 
representative of the total mm flux since the background 
averages to zero. The distribution of source sizes sup- 
ports the use of the 120" aperture, which is only smaller 
than 10 of the 55 BGPS sources. Figure [7] shows that for 
the 120" aperture, only approximately 5 source fluxes 
will be contaminated by emission from nearby sources. 

We have chosen to present the 5*^(120") fluxes as well- 
characterized fluxes for comparison with other studies. 
We present masses calculated from both 5*^(120") and 
5*1/ (int). However, when calculating physical properties 
which are highly dependent upon size, such as the volume 
and surface densities, we will present 5'!y(int) and use the 
derived source radius as the physical size. 

5.2. NH3 Results 

Example spectra of the NH3(1,1), (2,2), (3,3) and 
H2O lines are shown in Figure [8] for the source 
G189.776-I-00.343 (ID #13). Spectra for aU NH3 point- 



ings are shown in Figures 8b— at in the online materials. 
NIl3(4,4) and C2S were not observed for many sources 
and are not shown here. We detected a > 5cr NIl3(l,l) 
line in 95% (53 of 56) of all NH3 pointings, and 60% (33) 
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of the 46 cross-matched NH3 pointings and the peak of the asso- 
ciated BGPS source. The striped histogram marks the distance 
between the peak in 1.1 mm emission and the nearest NH3 point- 
ing for the 34 BGPS sources in the full sample. The NH3 point- 
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between pointings on order of the GBT beam size (31") or larger. 
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of the pointings had at least a Scr detection in both the 
(1,1) and (2,2) hnes, where a Scr detection means the in- 
tegrated intensity is greater than 5 times the measured 
error in the integrated intensity for each source. 

5.3. Combined 1.1 mm and NH3 Results 

Each NH3 pointing was assigned to a single distinct 
BGPS source. Of the 55 BPGS sources, 37 have ammo- 
nia observations along a line of sight that intersects the 
source. Three BGPS sources have an ammonia point- 
ing associated with the source that lacks any discernible 
NH3(1,1) emission. In the remaining 34 BGPS sources, 
there may be multiple NH3 pointings within the con- 
tours of the BGPS source, and the NH3 pointings do 
not necessarily correspond to either the centroid or peak 
coordinates of each BGPS source because they were cho- 
sen by eye prior to completion of the final BGPS maps. 
We exclude the seven NH3 pointings that do not fall 
within the contours of an extracted BGPS source from 
the analysis. These pointings correspond to low signifi- 
cance BGPS sources. The seven pointings all have signif- 
icant NH3(1,1) emission, indicating that many marginal 
detections in the BGPS not included in the catalog are 
likely real. 

In this analysis we will consider only the NH3 pointing 
which is located closest to the peak of the millimeter 
emission. The angular separation between the closest 
NH3 pointing and the peak of the millimeter emission 
is shown in Figure O The white histogram includes all 
NII3 pointings while the striped histogram includes only 
the NH3 pointings closest to the millimeter peak. Most of 
the NH3 pointings assigned to a BGPS source are located 
less than a FWHM beam away from the millimeter peak. 

We define a "full sample" which consists of 34 BGPS 
sources that have at least a 5a detection in the NH3(1,1) 
transition. Sources without a 5cr detection in the 
NIl3(2,2) transition can only provide an upper-limit to 



the derived gas properties such as Tr- . Consequently, the 
full sample can only characterize upper-limits for these 
properties. We also define "subset 1" which is comprised 
of 25 BGPS sources with a > 5cr detection in both the 
NH3(1,1) and (2,2) transitions. However, subset 1 may 
not be characteristic of the faint end of the full sample. 
Due to the complementary properties of each sample, we 
will present the characteristics of both. For clarity, we 
point out that subset 1 is entirely contained in the full 
sample. 

Observed properties of NH3 pointings in the full sam- 
ple are given in Table [5l including the source ID (column 
1), specific coordinates of the NH3 pointing which do 
not necessarily correspond to the peak of the millimeter 
emission (Columns 2 and 3), radial velocity and Gaus- 
sian 1/e line width (columns 4 and 5), and main beam 
temperature and integrated intensities for the NH3(1,1), 
(2,2), (3,3) and (4,4) transitions (columns 6-13). For un- 
detected lines, we present upper-limits for the main beam 
temperature and integrated intensities given as Tmb< 5 
RMS andW<5 RMS Av/^/iN), where RMS is the 
noise per channel of each spectrum, Av is the channel 
width in velocity, and N is the number of channels used 
to calculate the RMS. 

Table [6] lists derived gas properties including the mea- 
sured optical depth of the (1,1) transition (column 2), 
kinetic gas temperature (column 3; see §6.2p . excitation 
temperature of the NIl3(l,l) transition (column 4), de- 
rived sound speed and non-thermal velocity dispersion (a 
and (Tjvt; columns 5 and 6; see ^^, and flags for detec- 
tion or non-detection of the H2O and CCS lines (columns 
7 and 8). The additional 12 NII3 pointings are included 
in Tables [5] and |6l but are not included in the general 
analysis. Numbers given in parentheses are errors in the 
same units as each column. 

The distributions of the fiux densities, 5^(120") and 
5*1/ (int), are shown in Figure [TOl The white histogram 
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includes all 55 BGPS sources extracted from the mil- 
limeter maps. The gray histogram includes all BGPS 
sources included in the full sample (> 5a detection in 
NH3(1,1)), while the striped histogram denotes subset 
1 (> 5(7 detection in NH3(1,1) and (2,2)). The faintest 
sources are excluded from subset 1, as expected. 

Figure [11] plots the peak Tmb of the (2,2) transition 
versus the peak Tmb of the (1,1) transition. Pointings 
with < 5(7 detection in the (2,2) transition are shown 
as upper limits. The fainter lines in (1,1) are typically 
the pointings with a weak or non-detection of (2,2). The 
weaker BGPS sources typically have weaker NH3(1,1) 
detections and only upper-limits for the NH3(2,2) line. 

5.4. Maser Activity 

We also obtained observations of the 22 GHz water 
maser in conjunction with our GBT NH3 observations. 
Eight of the 34 BGPS sources presented here have wa- 
ter maser emission (see columns 7 and 8 of Table [6]) . Of 
these, only Sources #16 and #24 have not been reported 
previously as maser sources. The remaining 6 BGPS 
sources with maser activity were previously known to 
exhibit water as well as some SiO or methanol maser ac- 
tivity in a 1' region around the observed GBT position 
based on a SIMBAD search. The fraction, 8/34 = 0.2, 
can be regarded as a lower limit as two sources, #13 and 
#47, show maser emission along one line of sight but not 
another line of sight through the same core. Since the 
GBT ammonia spectroscopy was not complete mapping, 
some sources may have maser emission elsewhere in the 
object. However, the masing lines of sight are those near- 
est the peaks of the BGPS emission, so the lower limit is 
likely representative of the true fraction. For comparison, 
Szymczak et al. (2005) find a detection fraction of water 
masers of 0.52 using 6.7 Ghz methanol masers, associated 
with massive young stellar objects, as a prior. Similarly, 



Chambers et al. (2009) find detection rates of 0.16, 0.54 
and 0.59 for water masers in their IRDC cores classi- 
fied with Spitzer data as "quiescent," "red" and "active" 
respectively. The BGPS-selected water maser detection 
fraction is most similar to the detection fraction in "qui- 
escent" IRDCs determined by Chambers et al. (2009). 
We see no evidence for maser activity in the NH3(3,3) 
transition. All (3,3) spectra are consistent with the tem- 
peratures and line widths determined in conjunction with 
the (1,1) and (2,2) states. We defer further discussion 
of the water maser detection fraction to a later paper 
(M. K. Dunham et al., in prep). 

6. ANALYSIS 
6.1. Radial Velocities and Distances 

We adopt the VLBI parallax distance to S252 of 
2.10 ± 0.026 kpc (Reid et al. 2009) for the entire sam- 
ple of BGPS sources. The radial velocities observed in 
the NH3 spectra are listed in Table [5] (column 4), and 
are shown in Figure 1121 Filled circles mark sources in 
Subset 1, while the crosses mark sources without a 5a 
detection in the NH3(2,2) line. The sources excluded 
from Subset 1 are spread evenly in velocity space, as 
expected since the weaker BGPS sources should not be 
preferentially located at any given radial velocity. To- 
ward the Galactic anti-center (as well as all Cardinal 
points: I = 0°,90°, 180°, and 270°), radial velocities ex- 
pected from purely circular rotation are km s~^, ren- 
dering kinematic distances unsound. A non-zero radial 
velocity along these lines of sight is most likely due to 
the peculiar motions of individual gas clouds rather than 
galactic rotation. 

From Figure [T^ it is apparent that groups 2 and 3 
span the same velocity range (5 to 10 km s^^) while 
group 1 is found at slightly lower radial velocities (—1 
to 4 km s^^). Although there is no evidence in our data 
that these groups are related, surveys of lower density 
tracers (such as CO and ^^CO) toward S247 and S252 
(groups 1 and 2) show evidence of interaction between 
the two HII regions. Kompe et al. (1989) mapped the 
S247/252 complex in ^^CO ( J = 1 - 0) transition with a 
beam-sampled map at 4' resolution and found a bridge 
connecting the two clouds. Specifically, the bulk of the 
material towards S247 is found between velocities of — 1 
km s~^ and 5 km s~^ with a tail extending to veloci- 
ties of up to 10 km s^^, which overlaps the S252 region 
in velocity space (see Figure 7 in Kompe et al. 1989). 
Carpenter et al. (1995a) mapped a larger region includ- 
ing the S247 and S252 complexes in "C0(J =1-0) 
and found three distinct filaments. Two filaments with 
radial velocities of ~2 km s~^ wrap around S247, while 
the third filament, at a radial velocity of ^10 km s^^, 
extends to the S252 region. It is clear that the S247 and 
S252 regions (groups 1 and 2) are interacting, justifying 
our application of the VLBI parallax distance to our en- 
tire sample regardless of the difference in radial velocities 
between the groups of sources. 

Not all sources will be at precisely the parallax dis- 
tance; however, we account for variation in the distances 
by including a larger error than the parallax distance 
error. Since it is clear that groups 1 and 2 are part of fil- 
aments wrapping around and interacting between S247 
and S252, we will consider the size of the HII regions 
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when determining the spread in distance. S247 has a di- 
ameter of 9' and S252 has a diameter of 40' (Sharpless 
1959). If the cloud between Groups 1 and 2 is assumed 
to be spherically symmetric, its angular size on the sky 
can provide an estimate of the distance along the line of 
sight. Since it is apparent that groups 1 and 2 are inter- 
acting, we will consider a sphere of approximately 100' 
in diameter that encompasses both groups 1 and 2 and 
the associated HII regions, and corresponds to a distance 
uncertainty of 0.061 kpc. We will assume the spread in 
distance is approximately equal to the 100' diameter, and 
will assign that to the distance error to account for the 
case where the parallax source is found at the extreme 
edge of the region and note that this is likely a minimum 
estimate of the error. 

6.2. Kinetic Temperature 

The kinetic temperatures derived from the NH3 ob- 
servations are given in column 3 of Table |6l The un- 
certainties in individual determinations of Tk are small, 
with {'Jj'j. /Tk) = 0.04. In this and subsequent analy- 
ses, the uncertainties quoted are the standard deviation 
about the mean of the sample, rather than the uncer- 
tainty in the mean. For sources without a Scr detec- 
tion in NH3(2,2) the derived kinetic temperature is an 
upper-limit, and thus the mean kinetic temperature for 
the full sample is only an upper- limit. The distribution 
for the full sample is characterized by {Tk) < 19.4 ± 4.9 
K, while the distribution for Subset 1 is characterized by 
(Tk) = 20.6 ± 5.1 K. The mean limit on kinetic temper- 
ature for the sources without a 5ct NH3(2,2) detection is 
(Tk) < 16.1±2.3 K. As some property distributions may 
be non-Gaussian, the uncertainty should be regarded pri- 
marily as an indicator of the scatter in the sample. The 
dispersion about the mean kinetic temperature for the 
full sample is thus six times larger (0.24) than the frac- 
tional uncertainties in individual temperatures, so the 
differences among sources are real. 

Zinchenko et al. (1997) mapped 17 molecular clouds 
associated with FIR sources in NIl3(l,l) and (2,2) and 



found that Tk increases toward the peak of the NII3 
emission. Additionally, Friesen et al. (2009) have shown 
that NH3 emission closely follows mm continuum emis- 
sion. Since our NII3 pointings are located near the peak 
of the mm emission, which corresponds to the NII3 peak 
as well, the mean Tk is likely smaller than the peak Tk 
we report. 

6.3. Velocity Dispersion and Non-Thermal Pressures 

The line-width provides a measure of the one- 
dimensional velocity dispersion within each source. 
Rosolowsky et al. (2008) found velocity dispersions rang- 
ing from 0.07 to 0.7 km s"^ for the NH3(1,1) and (2,2) 
lines of low-mass star- forming cores in the Perseus molec- 
ular cloud, and they suggest that there is evidence of 
multiple velocity components in all sources with ay > 
0.2 km s^^. The velocity dispersions are listed in col- 
umn 5 of Table [5] and the distribution of our sample is 
shown in Figure [121 In contrast to the low-mass sources 
in Perseus, our distribution has a tail extending up to a 
high velocity dispersion a^ ~ 1.3 km s~^. 

The observed velocity dispersion is a combination of 
the thermal and non-thermal motions of the gas. In or- 
der to calculate the non-thermal contribution, a^T, we 
simply remove the thermal contribution: 



o'NT = \l crl 



kTK 
l7mH 



(4) 



where {UTk /^Trnn)^''^ is the thermal broadening due to 
NII3 and 17mH is the mass of a single ammonia molecule. 
The non-thermal velocity dispersion ranges from 0.37 
km s^^ to 1.3 km s~^ with a mean of 0.72 ±0.25 km s^^ 
for subset 1. We can compare (Jnt to the predicted ther- 
mal sound speed of the gas given by a = {kTK/f^rnHy/"^, 
where fi — 2.37. The thermal sound speed for subset 1 
ranges from 0.22 km s~^ to 0.34 km s~^ with a mean of 
0.27±0.03 km s~^. The mean Mach number is then 2.8. 
Both a and <7nt are listed in Table [6] (columns 5 and 6) . 
The corresponding ratio of ctatt to a ranges from 1.5 to 
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Fig. 13. — The non-thermal velocity dispersion, cr^fT, versus de- 
rived kinetic temperature. Arrows indicate upper limits for the 
kinetic temperatures. SoUd line denotes the thermal sound speed 
given by a = sJkT/ fimu . 

4.2 with a mean of 2.7±0.8. The non-thermal velocity 
dispersion is shown in Figure fT3l versus Tk- The solid 
line plots the thermal sound speed, a, as a function of 
Tk- In all sources, the non-thermal velocity dispersion 
is much larger than the thermal velocity dispersion and 
the discrepancy increases with Tk- 

There is no apparent correlation between (Jnt and 
source size or mass. However, there is a correlation be- 
tween the position of a source within its group and umt- 
Figure [14] displays the position of each source in Galactic 
coordinates and the symbol size denotes the magnitude 
of CTjvT- For Group 2, the highest values of ctatt are found 
at the center of the cluster. A similar statement could be 
made for Group 3, but there it is difficult to definitively 
say, due to the smaller number of sources. 

6.4. Virial and Isothermal Masses 

Combining heterodyne and dust continuum observa- 
tions provides estimates of both virial mass and total 
gas and dust mass. Total gas and dust mass (also called 
the isothermal mass, Mi^o, due to the simplifying as- 
sumption that the dust can be characterized by a single 
temperature) is given by 



13.1 Mp 







M,, 



Oi/ 



S,D^ 



K^B^{Td) 
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,13.0 A'/Td 



IJyJ \l kpcj \e^^-° ^/200 ^ - 1 



(5) 



where S^ is flux density, D is the distance, Ki, is the 
dust opacity per gram of gas and dust, and B^, is the 
Planck function evaluated at a dust temperature, Tj^. 
We logarithmically interpolate the Ossenkopf & Henning 
dust opacities (1994, Table 1, column 5, commonly re- 
ferred to as 0H5 dust) to the effective central frequency 
of the Bolocam bandpass convolved with a 20 K black- 
body modified by an opacity varying with frequency as 
K^ ex v^-^ (271.1 GHz; Aguirre et al. 2010) and obtain 



K„ = 0.0114 cm^ g-i. 

We assume the gas and dust are collisionally coupled, 
make the simplifying assumption that T^ = Tk and use 
Tk derived from the NH3 observations. Evans, Blair, & 
Beckwith (1977) observed the S255 molecular cloud in 
i^CO, i^CO, H2CO, and 2.2 /xm emission. They stud- 
ied the gas and dust energetics and concluded that the 
primary energy flow is through infrared dust emission 
and that an embedded source and the exciting stars of 
the many nearby HII regions are the likely heat sources 
for the cloud. Dust typically controls the energetics in 
the dense interior of the star-forming regions where the 
density is high enough to shield the gas from the ultra- 
violet radiation in the interstellar radiation field (ISRF) 
and also to transfer energy from the dust to the gas via 
colhsions (Goldreich & Kwan 1974, Evans et al. 2001). 
Towards the outer edges of the clump the density is much 
lower, and the gas is no longer collisionally coupled to the 
dust. The gas is also now subject to the ultraviolet por- 
tion of the ISRF, allowing the gas kinetic temperature to 
rise above the dust temperature (see Figure 9 in Mueller 
et al. 2002). The BGPS observations probe all of these 
regions, and detailed models will be needed to include 
temperature gradients. For now, we rely on previous 
modehng (e.g., Mueller et al. 2002, Shirley et al. 2002) 
to support the assumption that T^ — Tk on average. 

As discussed in ^16. 21 our single measurement of Tk 
likely corresponds to the maximum temperature seen 
within the BGPS source. Thus, using Td — Tk to de- 
scribe the dust temperature for the entire source will re- 
sult in a lower-limit for the isothermal mass. Zinchenko 
et al. (1997) find a maximum Tk of 28 K at the peak 
of the NH3 emission, and a minimum Tk of 15 K at the 
edge of the NH3 emission. This difference in tempera- 
tures corresponds to a factor of 2.3 increase in mass if 
the minimum temperature is assumed to be Td- If we as- 
sume mean values from Zinchenko et al. (1997) of Tr-=24 
and 17.5 K, the mass would increase by a factor of 1.5 as- 
suming the minimum temperature as Td rather than the 
maximum. If we are observing the maximum Tk then 
our calculated values of Miso are likely to underestimate 
the mass by a factor of up to approximately 2. 

The virial mass for a spherical, uniform density gas 
cloud is given by 



M„ 



5alR 



(6) 



where a^ is the Gaussian 1/e width of the line (model 
fit for the velocity dispersion), R is the source radius, 
and G is the gravitational constant. We assume that the 
velocity dispersion measured with our single NH3 point- 
ing is also applicable at the edge of the 1.1 mm con- 
tinuum emission, and use the single measured value in 
conjunction with the mm-derived object radius, Robj, to 
calculate the virial mass. This assumption is supported 
by Friesen et al. (2009) who found that NH3 emission 
closely traces 850 /im continuum emission and covers the 
same spatial extent as the cold dust producing the con- 
tinuum emission. However, Zinchenko et al. (1997) found 
that the ammonia line width increased toward the peak 
of the NH3 emission for about half of their sample. Thus, 
we are likely overestimating the virial mass by using the 
line width measured near the peak. 
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Myir, Miso{l'20"), and Afiso(int) are given in Columns 
2 — 4 of Table [7l For Subset 1, the isothermal masses 
within a 120" aperture range from 47 M© to 610 Mq 
with a mean of 280±180 M©. Similarly, Miso(int) ranges 
from 33 Mq to 850 Mq with a mean of 300±220 M©. 
The virial masses range from 33 Mq to 560 Mq with a 
mean of 220±160 Mq. The distributions of M,^o(120") 
and Miso (int) are shown in Figure [TSl The mean Alyir 
to Miso(int) ratio is 0.96±0.93 for Subset 1. The ratio of 
Myir to Miso(int) is shown versus Miso(int) in Figure [TCI 
The dashed line marks the mean ratio of 1.0. The trend 
of decreasing M^r / Miso reflects the dependence of the 
ratio on the inverse of the source size at a given distance. 
Robj increases with Afiso(int) which lowers Myir / Miso as 
a function of Miso(int). 

In these calculations we have made the simplifying as- 
sumption that the density is uniform. However, if the 
emission is not spherical or if the density follows a power- 
law the virial mass will differ. A corrected virial mass is 
given by 
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for p < 2.5, 



(7) 
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(Bertoldi & McKee 1992) where ai is the correction for a 
power-law density distribution, 02 is the correction for a 
non-spherical shape, and p is the density power-law index 
from n{r) = nf{r/rf)-P (Mueller et al. 2002). The mass 
ratio is then given by 
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The mean aspect ratio of our sources is {Rmaj / Rmin) — 
1.5 ± 0.3, which gives 02 = 1.1. To estimate the ef- 
fect of a power-law density distribution, we adopt the 
mean power-law index found by Mueller et al. (2002), 
{p) — 1.8 ± 0.4, which was derived by modeling the 
emission from 31 massive star-forming regions, and find 
ai = 1.43 ± 0.34. If the corrections for a power- law den- 
sity distribution (Equation [S]) and for a non-spherical 
shape are included, then {Mmr/Miso) would be decreased 



by 36% to 0.61±0.59 for Subset 1. However, since we 
do not have accurate values of the power-law index for 
our individual sources, we will not include the correction 
in our calculations or subsequent discussion. Magnetic 
fields are also excluded from our virial mass calculation, 
but if magnetic fields provide additional support to the 
cloud, then the true mass would be larger than given by 
Equations |6] or [T] 

The observed scatter in Myir to Miso (Figure [T6|) could 
be a result of deviations from the assumed distance of 
2.10 kpc or deviations from Equation IH] due to differ- 
ent power-law density distributions or source ellipticities, 
all of which could differ from source to source and pro- 
duce the observed scatter. Additional assumptions which 
could produce the observed scatter include the assump- 
tions that the sources are virialized and that the velocity 
dispersion measured at the peak of the mm emission is 
also valid at Robj- A virial parameter, Myir/ Miso, of 1 
suggests that both mass estimates are reasonable. 

6.5. Volume and Surface Densities 

Average particle density and average surface density 
can be calculated based on a mean particle weight per 
H nucleus (/J,p=2.37) and a total mass within a given 
physical size, assuming a uniform density. We assumed 
a spherical volume given by 47ri?j^j, /3. Mean number 
densities of particles and mass surface densities were cal- 
culated for each source (see Column 5 of Table [71) . For 
Subset 1, np(int) ranges from 1.4 x 10^ cm~'^ to 17.0 x 10^ 
cm~^ with a mean of 6.2 ±4.4 x 10^ cm~'^, and the char- 
acteristic values for the full sample are given in Table 

e 

The distribution of derived surface densities, I](int), 
ranges from 0.019 g cm^^ to 0.13 g cm^^ with a mean 
of 0.055±0.032 g cm-^ for Subset 1. The mean sur- 
face density for the full sample is 0.048±0.030 g cm^^ or 
230±140 Mq pc~^ which is similar to the gas surface den- 
sity threshold for star-formation of 200 Mq pc~^ seen in 
nearby, low-mass star-forming clouds (A. Heiderman et 
al., in prep.). This similarity supports the idea that our 
BGPS sources are sites of star formation, while the low 
volume densities imply we are detecting the less dense, 
larger scale portions of clumps rather than the higher 
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Fig. 15. — Distributions of derived Migo derived from 5^(120") (left) and S'^(int)(right). The white histogram marks the full sample and 
the striped histogram marks subset 1. 



density cores. 

There is also evidence that a few of the BGPS clumps 
in this region include multiple higher density sources 
when observed with interferometers (see SJT]). Eight of 
the 34 BGPS cores also have water maser emission (see 
Tableini) which requires densities greater than lO^'^ cm~'^ 
(Strelnitskij et al. 1984; Elitzur et al. 1989), suggesting 
that the clumps contain much higher density substruc- 
ture. 

6.6. Excitation Temperature 

The measured excitation temperatures are plotted 
against the kinetic temperatures in Figure [T71 The filled 
circles represent Subset 1, and the corresponding distri- 
bution is characterized by (Tcx(l,l)) = 7.8±5.6 K. There 
are two sources where the best fitting model requires the 
limiting case Tc^=Tk, and higher quality data would be 
needed to constrain Tex- 

The mean difference between kinetic and excitation 
temperatures for Subset 1 is {Tk-T^^) = 12.8 ± 4.5 K. 
There are two possible explanations for the differences in 
derived temperatures. The first possibility is inherent in 
the way we measure the excitation temperature. As de- 
scribed in §4.21 in order to calculate Tox we must assume 
something about the ratio of the solid angle of the source 
to that of the beam, which here we set to unity. A lower 
excitation temperature could be artificially induced by 
over-estimating the beam filling factor. 

Another possible explanation is the density in the re- 
gion. The level populations are determined by both 
coUisional and radiative processes which depend on the 
local density. Pavlyuchenkov et al. (2008) study the 
conditions under which molecular rotational lines form 
using radiative transfer simulations of C0(J=2— 1) and 
HC0+(J=1— 0), and define a thermalization density nth, 
where Tex is within 5% of Tr-. At densities greater than 
nth, collisions regulate the level populations by both ex- 
citing and de-exciting molecules. A commonly used mea- 
sure is the critical density, defined by ncrit = ^/7, where 
A is the spontaneous emission rate and 7 is the colli- 



sion rate coefficient. For low frequency transitions, such 
as the NH3 inversion lines, collisions are negligible for 
n < ncrit, and the excitation temperature is set by the 
cosmic background radiation. At intermediate densities 
where ncrit < n < nth, the excitation temperature is 
set by both collisions and radiation and is equal to nei- 
ther the background radiation temperature nor the ki- 
netic temperature. 

Using an LVG code (Snell et al. 1984), we calculate Tgx 
as a function of density for a gas temperature of Tk= 20 
K, with A^(NH3) = 10^^'^ cm~^, over a density range of 
100—10^ cm""^. We assume a line width of 0.75 km s~^ 
to include radiative trapping. Figure [TH] plots Tex as a 
function of volume density for the NH3(1,1) transition. 

From these models we find nth — 210ricrit for the 
NH3(1,1) transition. The large difference in Ucrit and 
nth for the NH3 transitions is due to the low Einstein 
A coefficient and the importance of stimulated processes 
to the level populations. For comparison, Evans (1989) 
finds n{Tex= 0-9Tk) = l■33ncr^t for the CS(3-2) tran- 
sition. Although this criterion is less stringent than 
nth, it demonstrates the differences between the low fre- 
quency radio transitions and higher frequency submil- 
limeter transitions. The low frequency transitions, such 
as the NH3 transitions, require a density much higher 
than the critical density before Tox begins to approach 
Tk- 

Although the volume-averaged densities derived from 
the millimeter continuum are approximately equal to the 
critical densities for both the NH3(1,1) and (2,2) transi- 
tions (1.8 X lO'^cm"'^ and 2.1 x 10^ cm~^ respectively; 
Tielens 2005), the densities are not high enough for the 
gas to be completely thermalized, and Tcx<Tk as ex- 
pected at densities near the critical density. 

6.7. NH3 Excitation Density 

If we assume the gas is not in LTE, we can define the 
excitation density, Ugx, as the density required to pro- 
duce level populations described by Tex for a given gas 
kinetic temperature. We calculate Uex from Tk, Tex, and 
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Fig. 16. — Virial parameter {M^ir/ Miao) versus A/iso(int). The 
dashed Une marks the mean virial parameter of 1. 



r via the balance of a two-level system. The excitation 
density is given by (Foster et al. 2009; Swift et al. 2005; 
Casein et al. 2002) 



where 



k{J{Tcy.) ~ JjTcrnb)) 
/li.(l,l)(l-J(Tcx)/J(rK)) 



ricritP, (10) 



JiT)^ 



/c(l-e-'''^(i'i)/'=^)' 
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Tcmb = 2.73 and the escape probability is given by /? = 
(1 — e~'^)/T, where we take r = T(i^i) x 0.233 which is 
the maximum optical depth in a single hyperfine line of 
the NH3(1,1) transition. 

The excitation densities are listed in Column 6 of Table 
[71 We cannot calculate n^^ for the two sources where we 
set Tex=TK (see i 36.6|) . The excitation densities for Sub- 
set 1 range from 2.1 x 10^ cm~^ to 29.9 x lO'^ cm~^ with 
a mean of 9.6 ± 7.1 x 10'^ cm"'^. Figure [T9l shows a com- 
parison of the NH3 excitation density and the volume- 
averaged density (rip) from the millimeter data. The 
solid line marks a one-to-one correlation, and is not a 
fit to the data. Although there is considerable spread, 
Hp is in good agreement with the excitation density. The 
mean ratio of np(int) to excitation density is 1.7 ± 1.0 
with a median of 1.8 for subset 1. The agreement be- 
tween Hex and Hp provide additional support that the 
NH3 and continuum emission trace the same physical re- 
gion. 

6.8. Column Densities 

The column density per beam of NH3 in the (1,1) state 
can be calculated from the observed parameters via the 
following equation: 



7V(1,1) = 2 



(1,1) 51 



1 



52 ^(14 



T{v)dv, 
(12) 




Fig. 17. — Excitation temperature of the NH3(1,1) transition 
versus kinetic temperature for the full sample. The filled circles 
mark subset 1, and crosses mark sources with upper limits in Tj^jj^. 
The dashed line marks a one to one correlation but is not a fit to 
the data. 

where 



T{v)dv —T [1,1) I dv y Si exp 
=V27r— J^(i.i)t(1,1) 



(13) 



and 1^(1.1) is the central frequency of the (1,1) transi- 
tion, T(i 1) is the total optical depth in the (1,1) transi- 
tion, gi and 172 are the statistical weights of the upper 
and lower levels of the (1,1) inversion transition and are 
equal, ^(1,1) = 1-68 x 10~^ s~^ is the Einstein A coeffi- 
cient for the transition (Pickett et al. 1998), and a^ is 
the Gaussian 1/e line width. The leading factor of 2 in 
equation [T^] assumes an equal number of NH3 molecules 
in each state of the inversion transition, which is justi- 
fied since the states have equal statistical weight and are 
separated in energy by 0.05 K. The total optical depth in 
the (1,1) transition is split over the 18 hyperfine compo- 
nents, which is expressed in the summation in Equation 
1131 Si is the weight of each hyperfine component, Vi is 
the frequency of each hyperfine component, and Oi is the 
Gaussian 1/e width of each hyperfine component. 

The total NH3 column density, Nmh^ can be calcu- 
lated by scaling A^(l, 1) by Z/Z(J = 1) where Z is the 
partition function given in Equations [5] and |31 Z{J = 1) 
is the population in only the J = 1 level, and we truncate 
the summation at 50 terms. The NH3 column densities 
are given in column 8 of Table [71 Nmh3 ranges from 
4.2 X 10" cm~2 to 75.8 x 10" cm"^ with a mean of 
3.5 X 10^"* cm~^ for subset 1. 

We can also calculate the column density of H2 via 
the millimeter dust observations. From the gas and dust 
mass within a given aperture, we can calculate Nh2 by 



Nh,^ 



-^*J-iso 



fj.H^rriHTrRlhj' 



(14) 



where ij,h2 — 2.8. Nh2 is listed in column 9 of Table 
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Fig. 18. — Excitation temperature versus volume density for 
NH3(1,1) for Tk=20 K. Vertical dashed line denotes n = 
"crit(l,l)i and the dotted line denotes the mean particle density 
found for subset 1. The LVG code predicts an excitation tem- 
perature of 6.3 K at the mean density of our sample. Our mean 
excitation temperature is 7.3 K. 

[3 Nh^ ranges from 4.1 x 10^^ cm-^ to 27.1 x 10^^ 
cni~^, with a mean of 1.2 ± 0.7 x 10^^ cm~^. Using 
the relationship between H2 column density and visual 
extinction of Bohhn et al. (1978) assuming Ry = 3.1, 
Nh^/Av = 9.4 X 10^° cm-2 mag~\ we find Ay = 12 A 
mag for the mean H2 column density of our sample. How- 
ever, Rv — 3.1 may not be an accurate description of 
the dust in these dense regions. Chapman et al. (2009) 
find that regions with Ak > 1 mag {Ay > 5 mag) are 
described well by a Ry = 5.5 extinction model (Wein- 
gartner & Draine 2001). Using this extinction model, 
NH2/AV = (2 X 1.086Cext)~^: where the factor of two 
converts Nh to Nh2 and Cext is the extinction cross sec- 
tion in the V band. Using Nh^/Av = 6.86 x 10^° cm'^ 
mag^^ we find Ay — 17.5 for the mean column density 
of our sample. 

We can obtain a measurement of the NH3 abundance 
by comparing the column densities of NH3 and H2 . The 
abundances are presented in column 10 of Table [T] The 
abundance ranges from 0.89 x 10~* to 8.3 x 10~^ with 
a mean of 3.1 x 10~* for subset 1. These results are 
similar to the values derived for dense cores, e.g. Harju 
et al. (1993, 3 x 10^^), Tafalla et al. (2006, 2.8 x 10"*) 
and Foster et al. (2009, 2 x 10"^). 

7. DISCUSSION 

The physical properties of the BGPS sources in the 
Gemini OBI Molecular Cloud show that the majority 
of the sources are the clumps from which entire stellar 
clusters will form. Characteristics of all properties are 
listed in Table H] for both the full sample and subset 1 , 
where numbers in parentheses are the standard deviation 
rather than an error in the mean. Recall that clumps 
are typically characterized by masses ranging from 50 -— 
500 Mq, radii of 0.3 — 3 pc, mean densities of 10'^ — 
10^ cm~^, and gas temperatures of 10 — 20 K (Bergin 
& TafaUa 2007). In the Gemini OBI region, we find a 
mean isothermal mass of (Miso(int)) = 200 ± 220 M0, 
a mean radius of (R) = 0.56 ± 0.19 pc, a mean density 
of (np(int)) = 6.2 ± 4.4 x 10'^ cm~'^, and a mean gas 
kinetic tempearture of (Tk) = 20.6 ± 5.1 K. The mean 
properties of our BGPS sources fall squarely in the ranges 
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Fig. 19. — NH3 excitation density, riex, versus np{int). The 
dashed line marks a one-to-one correlation and is not a fit to the 
data. The mean volume densities are in good agreement with the 
NH3 excitation densities. 



of characteristic clump properties. 

In contrast, the unresolved BPGS source in this re- 
gion (ID #51) has properties within the range of the 
characteristic properties of cores. Recall that cores have 
masses from 0.5 — 5 M©, sizes of 0.03 — 0.2, mean den- 
sities of 10* — 10^, and gas temperatures between 8 
and 12 K (Bergin & Tafaha 2007). Source #51 has 
Afiso(int)= 44±13 Mq, an upper-limit for size of R< 0.17 
pc, a mean density of np(int)= 2.9 x 10"^ cm~'^, and an 
upper-limit of gas kinetic temperature of Tk< 14.8 K. 
Afiso(iut) is significantly above the range of core masses, 
but also falls below the range of characteristic clump 
masses. The mean density does not fall within the char- 
acteristic density range of cores, but the upper-limit used 
for the source size could be the explanation. An unre- 
solved source must have a size smaller than the beam 
size, so the beam size is artifically giving a lower den- 
sity. There are additional unresolved BPGS sources but 
they do not have corresponding NII3 pointings and were 
excluded from our analysis. 

Further evidence that the BGPS sources are not di- 
rectly forming a single massive star is in the mean sur- 
face density of the sample, {Yj{int)) = 0.055 ± 0.032 g 
cm~^. Krumholz & McKee (2008) have shown that a 
surface density of at least 1 g cm~^ is required to avoid 
fragmentation and form a single massive star. Although 
the mean masses of our clumps are high enough to form 
a massive star given typical star formation efficiencies, 
the mean surface density is not high enough to prevent 
fragmentation and suggests that there should be signifi- 
cant substructure within these clumps. Substructure has 
been seen in higher resolution studies of some of the Gem- 
ini OBI BGPS clumps and is discussed in detail in the 
following section. Additionally, the mean virial parame- 
ter for Subset 1 is approximately 1 which suggests that 
the BGPS sources are stable against collapse on the size 
scales the survey is sensitive to. 
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7.1. Comparison to Other Studies 

A comparison to the literature provides further evi- 
dence that the BGPS sources are hkely the clumps from 
which stellar clusters will form. The BGPS sources often 
resolve into multiple smaller, higher density regions when 
observed at higher resolution with telescopes such as the 
SMA, SCUBA on the JCMT, and SHARC on the CSO. 
Densities calculated from 1.3 mm continuum SMA obser- 
vations are on the order of 10® — 10^ cni^^ (Cyganowski 
et al. 2007), a factor of 10"^ greater than the densities 
of the BGPS sources presented here. Models of observed 
spectral energy distributions require both a hot (^ lOOK) 
and warm (~ 40K) dust component (Minier et al. 2005), 
again suggesting that the BGPS sources are the cooler, 
lower density clumps which contain hotter cores. When 
comparing our 1.1 mm emission to the large-scale CO 
emission, we find that the BGPS sources have H2 col- 
umn densities a factor of 10 greater than the mean col- 
umn density of the regions of the molecular cloud without 
star formation (Carpenter et al. 1995a). The resolved 
BGPS sources are typically clumps which are embedded 
in a molecular cloud and contain higher density, warmer 
cores while the unresolved BGPS sources are likely to be 
cores. Detailed comparisons for a few sources and the 
entire molecular cloud are presented below. 

7.1.1. S252 

S252 corresponds to G188.948-H00.883 (ID #26) in the 
BGPS. Minier et al. (2005; hereafter M05) observed 

methanol and CO toward S252, and also modeled mid- 
IR to millimeter wavelength emission. Using SCUBA on 
the James Clerk Maxwell Telescope and SIMBA on the 
Sweedish ESO Submillimetre Telescope (SEST), they re- 
solve S252 into two millimeter sources, G188.95-MM1 
and — MM2. Models fit to the spectral energy distri- 
butions require a warm extended envelope (42 and 50 
K) for each source and a hot central component (150 
K) for MMl. The rotational temperature derived from 
their methanol observations is 58 K. The temperatures 
derived from both SED fitting and molecular line obser- 
vations are all significantly warmer than our NH3 kinetic 
temperature of 29 K. In addition to warmer tempera- 
tures, M05 also derived a smaller mass than the BGPS. 
The BGPS does not resolve two sources here, and finds 
Miso (int) =556 M©, while M05 finds a total mass for 
MMl and MM2 of 155 M© when scaled to our assumed 
distance of 2.1 kpc. 

7.1.2. S252A 

Mueller et al. (2002; hereafter M02) observed S252A 
in 350/im continuum emission using SHARC at the CSO, 
and modeled the observed SED including wavelengths 
ranging from 12 /im to 1.1 mm. They fit a power-law 
density distribution and determined a density normalized 
at a radius of 1000 AU, n/. S252A in the M02 sample 
corresponds to G189.776-H00.343 (ID #13) in our sam- 
ple. M02 calculated an isothermal dust mass of 220 M© 
within a 120" aperture when scaled to our assumed dis- 
tance. We find Mj^o(120")= 373 ± 40Mofrom the 1.1 
mm emission. SHARC filtered out extended emission on 
a smaller scale than the BGPS, so the M02 observations 
detected the warmer, higher density, and lower mass sub- 
structure within the BGPS source. The Tiso given by 



M02 was 30K, marginally warmer than the temperature 
derived from the NH3 observations, Tk— 26K. 

In general, the M02 study included higher-mass, 
warmer sources than our sample. Their sample selec- 
tion required the presence of strong water maser emis- 
sion (Arcerti water maser catalog; Cesaroni et al. 1988) 
as a signpost of massive star formation. M02 found 
(M,^o(120")) = 2040 ± 4410 Mq, with a median of 400 
Mq for their sample of 51 massive star-forming regions. 
In contrast, we find (M,so(120")) = 280 ± 180 M©, with 
a median of 190 M©. Our sources are approximately a 
factor of 10 less massive than the sample studied by M02. 

7.1.3. S255N and S255IR 

S255N corresponds to G192. 581-00.043 (ID#45) in our 
sample, while S255IR corresponds to G192.596-00.051 
(ID #47). We find M,^o(int)= 397 ± 32 M© and Tk= 
30K for S255N, and Mj^o(int)= 540 ± 43 and Tk^ 32K 
for S255IR. We compare our BGPS and NH3 derived 
properties with studies by Zinchenko et al. (1997; here- 
after Z97), M05, and Cyganowski et al. (2007; hereafter 
C07) in detail. 

Z97 mapped NH3(1,1) and (2,2) emission from S255 
with the 100 m Effelsberg radio telescope. They cen- 
tered their map on an IRAS point source located be- 
tween S255N and S255IR, and considered the emission 
from both to be part of a single source. They report 
observed values at two locations, one of which corre- 
sponds very well with S255N and our source G192.583- 
0.043. The observed NH3 properties between the two 
studies agree fairly well considering the NH3 positions 
differ slightly. For example, Z97 find visr — 8.65 km s~^ 
and a^ = Aw/(8 In 2)^/2 = 1.02 km s'^, while we find 
visr — 8.66 km s~^ and a^ = 1.31 km s~^. Z97 find that 
in half of their sources the NH3 line width increases to- 
ward the NH3 peak which could account for the same 
radial velocity but different line widths. Both studies 
derive slightly different Tk, with Z97 finding Tk^23 K 
and our study finding Tk^SO K for S255N and 32 K for 
S255IR. The Z97 mass derived from the NH3 column den- 
sity and a canonical abundance is 420 M© when scaled 
to our assumed distance. Since this includes both S255N 
and S255IR, this mass should be compared to the sum 
of masses for the two BGPS sources, Miso(int,tot) = 937 
M©. Our total M^so (int) is comparable to the virial mass 
Z97 derive, Myir{Z97) = 1134 M© when scaled to the 
maser parallax distance. 

C07 observed S255N in 350 /im continuum with 
SHARC on the CSO and 1.3 millimeter continuum with 
the SMA. S255N-MM1 corresponds to G192. 581-0.043 
in our sample (ID #45). SHARC {Q^h =15") recov- 
ers only a single peak at the position of S255N-MM1, 
while the SMA (4.7" x 2.4" beam) clearly resolves three 
separate sources: SMAl, SMA2, and SMA3. The total 
beam-averaged mass of all three sources scaled to the 
parallax distance is 17—42 M©, where the range is due 
to uncertainty in the assumed dust temperature. The as- 
sumed temperature range for SMAl is 40—100 K, while 
the other two sources are likely between 20 and 40 K 
based on the presence of molecular emission in SMAl and 
the lack thereof in SMA2 and SMA3. We find Tk= 29.5 
K for source #45 in our sample, which is consistent with 
the range of temperatures for SMA2 and 3 and is just 
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below the lower end of the temperature range for SMAl. 
The difference in the temperatures reflects the different 
material probed with the two sets of observations: the 
SMA observes the warm, high density structure within 
the cooler, lower density clumps to which the BGPS and 
GBT observations are sensitive. 

C07 derived densities in the range of 10^ — 10^ cm~'^ 
from the SMA 1.3 mm continuum emission. These den- 
sities are ^ 10'^ times higher than the derived BGPS 
densities. With a smaller beam size, the SMA resolves 
out the low density material which the BPGS is sensi- 
tive to and is able to detect the smaller, higher density 
regions within the BGPS clump. Our calculated Up is 
a volume averaged density which will underestimate the 
peak density of the regions where massive stars will form 
because a larger percentage of the volume is at the lower 
densities. 

Similar to C07, M05 found that S255N requires a hot 
component of 106 K and a warm component of 44 K 
while S255IR requires a hot component of 225 K and 
a warm component of 42 K based on fitting models to 
the observed spectral energy distribution. M05 find Miso 
of 140 and 210 Mq for S255N and S255IR, respectively. 
Their masses are significantly smaller than our isother- 
mal masses, which is consistent with the paradigm in 
which observations with a smaller beam size will resolve 
out the lower density emission and detect the substruc- 
ture within the BGPS sources. 

Chavarria et al. (2008) provides direct evidence that 
S255IR and S255N are forming a cluster. They used 
Spitzer IRAC and NIR data to identify sources with 
near- or mid-IR excesses. They identified a total of 10 
clusters of YSOs as well as a population of distributed 
YSOs. The largest cluster they identified corresponds to 
both S255IR and S255N (ID#45 and #47 in our sam- 
ple). They detected 140 stars with an IR excess, and 
assume the total cluster membership to be 488 sources 
after accounting for completeness and sources with no IR 
excess. 

7.1.4. S258 and G192.719+00.043 

There is also direct evidence of clusters forming in both 
S258 (G192.644-H00.003; ID #48) and G192. 719+00.043 
(ID #49) Chavarria et al. (2008) identified 49 sources 
with IR excesses in S258, and estimate the total num- 
ber of cluster members to be 186 after considering com- 
pleteness limits and the number of potential sources 
without an IR excess. The cluster they identified in 
G192.719-f00.043 includes 22 detected members and is 
estimated to include 44 total members. 

7.1.5. Gemini OBI Molecular Cloud 

Carpenter et al. (1995a; hereafter C95a) surveyed 32 
square degrees of the Gemini OBI molecular cloud com- 
plex in 12C0(J = 1-0) and ^^CO{J = 1-0). The 
detailed CO maps yield a wealth of physical properties, 
which provide an independent check of the properties 
determined in this study. C95a separated the molec- 
ular cloud into seven subregions marked by bright ex- 
tended ^^CO emission. Subregions 1—3 are associated 
with optical HII regions and correspond to groups 1-3 in 
this paper. Subregion 6 contains a Herbig Ae/Be star 
HD250550, but is not included in our BGPS image. The 



remaining subregions have no optical evidence of star 
formation. 

Along the lines of sight with a 5a detection in ^'^CO, 
C95a assume the ^^CO is optically thick and can deter- 
mine Tk- For subregions 4—7 with no star formation, Tr- 
is typically 6 — 8 K. In contrast, in subregions 1 — 3 with 
active star formation and optical HII regions the bulk of 
the material is found to be 6 — 8 K but the distribution 
has a warm tail reaching up to Tk= 25 K. The kinetic 
temperatures (~ 10 — 30 K) we derive from the NII3 
spectra are consistent with the warmer CO lines of sight 
where active star formation is likely occurring. C95a also 
determines the column density of II2 along lines of sight 
with a 5a detection of ^^CO. The subregions without star 
formation have {NH2) = 1.0 — 1.3 x 10^^ cm~^, while the 
cloud as a whole ranges from 1.2 x 10^° cm^-^ to 6.4 x 10^^ 
cni^'^. This is again consistent with our study, where we 
find {Nn^iint)) = 1.2 x 10^^ cm'^. 

For the subregions which correspond to our groups, we 
can determine the fraction of mass in dense clumps. To 
do so, we compare Mi^o (int) with Mlte from C95a. The 
total mass of dense millimeter clumps is 1835, 4549, and 
1973 M0 for groups 1 — 3, respectively. The correspond- 
ing Mlte for subregions 1 — 3 are 3.5 x 10'*, 5.2 x 10^, and 
4.3 X 10" Mq. We find 5.2%, 8.7%, and 4.6% of the total 
mass is in the dense clumps for subregions and groups 
1 — 3, respectively. C95a find that 10% of the molecular 
mass traced by CO is contained in lines of sight with II2 
column densities greater than 10^^ cm~^, similar to the 
mean column density of our clumps. This percentage is 
also similar to nearby, low-mass, star-forming molecular 
clouds Perseus, Ophiuchus, and Serpens, where < 10% 
of the total mass is in dense cores (Enoch et al. 2007). 
The fraction of mass found in the dense regions asso- 
ciated with star formation is very similar using different 
tracers as well as in high-mass and low-mass star-forming 
regions, two very different environments. 



7.2. What 



BGPS 



The BGPS is detecting thousands of sources which 
span a vast range of distances, and it is important to un- 
derstand the biases and trends inherent in the data. In 
contrast to the inner Galaxy BGPS sources, the Gemini 
OBI molecular cloud provides an excellent opportunity 
to characterize the clump properties without large un- 
certainties in distances. Two neighboring sources in the 
inner Galaxy could be many kpc apart along the line of 
sight, whereas the Gemini OBI sources are all located at 
the same distance. Since we have an accurate measure 
of the distance to the Gemini OBI sample we will use it 
as a starting point to understand the effects of distance 
on the derived properties. 

The detection criterion for a BGPS source is equiva- 
lent to a 5a detection within an area the size of a beam 
(Rosolowsky et al., 2009). This results in a minimum 
mass per beam, below which the emission will be lost 
in the noise or removed during the sky cleaning pro- 
cess. This minimum mass is given by Equation [SI with 
Si,=5a, where a is the RMS noise per pixel in units of 
Jy bearn"^ in the BGPS image, D is the distance, k^ is 
the dust opacity per gram of gas and dust, and B^{Td) is 
the Planck function evaluated at the dust temperature. 
Thus the BGPS is subject to the Malmquist bias: as 
the distance increases, the minimum mass to which the 
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BGPS is sensitive also increases. 

The minimum volume-averaged particle density to 
which the BGPS is sensitive is given by 



15cr 



1 



''P,rmn 



ATTfimHK,B,{T,)Ql,D el,D' 



(15) 



where fi = 2.37, mn is the mass of a hydrogen atom, and 
&mb is the FWHM of the CSO beam at 1.1 mm (33"). 
As the distance increases, the minimum mass increases 
(Equation[5]), and the minimum density decreases (Equa- 
tion lisp for sources smaller than the beam size, quanti- 
fying how BGPS observations are subject to Malmquist 
bias and beam dilution respectively. Similarly, if we allow 
the aperture size for photometry to vary at a constant 
distance, larger sources will have lower densities. 

For the Gemini OBI molecular cloud observations, a ~ 
0.37 Jy/beam, requiring S^(heam.)= 0.35 Jy for a 5a 
detection. The corresponding minimum mass is Mmin = 
20 Mq at a distance of 2.1 kpc and rD=20 K. Thus the 
minimum density required to detect a beam-sized source 
after deconvolution in the Gemini OBI molecular cloud 
is np^min ~ lO** cm~^. Sources larger than the beam size 
can have densities lower than 10'' cm""^. The smallest 
BGPS sources could have the highest densities and may 
be the cores which will form individual stars, multiple 
star systems, or small groups, whereas the larger sources 
will have generally lower densities and are likely to be 
the clumps which will form stellar clusters. 

In the inner Galaxy, the largest distance will greatly 
exceed that of the Gemini OBI molecular cloud and the 
minimum density will be significantly lower. At distances 
greater than approximately 2 kpc, the BGPS is only able 
to detect clumps which will form stellar clusters, and 
not individual cores. Since the Gemini OBI molecular 
cloud lies at 2.1 kpc, we may detect some cores as well 
as clumps. 

8. SUMMARY 

We present the image and catalog of the Gemini OBI 
Molecular cloud complex from the 1.1 mm Bolocam 
Galactic Plane Survey. We observed NH3(1,1), (2,2), and 
(3,3) inversion transitions and 22 GHz H2O maser tran- 
sitions toward 37 BGPS sources, as well as NH3(4,4) and 
C2S toward 10 sources with the GBT. We have obtained 
a detection > 5a in the (1,1) transition of 34 of the 37 
BGPS sources, and a detection > 5a in the (1,1) and 
(2,2) transitions in 25 of the 34 BGPS sources. 

We use a seeded watershed algorithm to extract the 
millimeter continuum sources from the BGPS maps. We 
find a mean FWHM of sources of 110", or 1.12 pc at a 
distance of 2.1 kpc. We estimate the gas physical prop- 
erties by fitting a simple model to all NH3 transitions 
simultaneously. From these fits, we find (Tk) = 20 K, 
and velocity dispersions ranging from 0.4 — 1.3 km s^^. 
We calculate the non-thermal velocity dispersion, com- 
pare it to the expected thermal sound speed, and find a 
ratio of cr^vT to a ranging from 1.5 to 4.2 with a mean 



of 2.6 ± 0.8. The velocity dispersions of all sources are 
dominated by non-thermal motions. 

We calculate the isothermal mass for each millimeter 
source based on S'^(120") and 5,y(int) assuming T^ = Tk 
from the NH3 observations. We find (Mj5o(120")) = 
230 ± 180 M0 and (Af,so(int)) = 250 ± 210 Mq, where 
the error is the standard deviation of the sample. We 
also calculate the virial mass within the source size de- 
rived by the source extration algorithm under the as- 
sumption that the NH3 emission traces the millimeter 
emission weU. We find {My^r) = 190 ± 160 Mq. The 
mean Myi^ to Miso ratio is 1.0±0.95. The total isother- 
mal mass for the portion of the Gemini OBI molecular 
cloud we have mapped is 8,400 M©, which represents 
6.5% of the total mass seen in CO observations. 

We also calculate the volume-averaged densities and 
find (np(int)) = 5.9 ± 4.3 x 10^ cm-^^ with the high- 
est density being np^max{^^t)= 1-7 x lO** cm~^. The 
mean surface density of the full sample is (E(int)) = 
0.048 ± 0.030 g cm-2. We calculate both NH3 and H2 
column densities and find mean values of 3.0 x lO^'' cm~^ 
and 1.0 X 10^^ cm~^, respectively. We find a mean NH3 
abundance of 3.0 x 10~^ relative to H2. 

The mean sizes, densities, kinetic temperatures and 
surface densities of the BGPS sources suggest they are 
cluster-forming clumps that contain the cores that will 
form individual or small multiple systems. 
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TABLE 5 
Observed NH3 Properties 



ID 

Number'' 



RA 
(J2000) 



Dec 

(J2000) 



(km s^'-) 



a(l,l_) 
(km s ^) 



Tmb(i.i) 

(K) 



W(l,l) 
(K km s-i) 



Tmb(2>2)'^ 
(K) 



W(2,2)b=<; 

(K km s-l) 



Tmb(3,3)'''= 
(K) 



W(3,3)'^''^ 
(K km s-i) 



Tmb(4,4)'^-= 
(K) 



W(4,4)b.c 
(K km s-i) 



01 06 07 30.8 

03 06 07 47.7 

04 06 08 18.2 

05 06 08 19.8 

06 06 08 23.5 

07 06 08 25.3 

11 06 08 33.3 

12 06 08 31.0 

13 06 08 35.2 

B 06 08 29.3 

16 06 08 40.6 

17 06 08 40.4 

18 06 08 42.1 

19 06 08 43.4 

20 06 08 43.9 

21 06 08 44.4 

B 06 08 47.0 

22 06 08 45.9 

24 06 08 49.5 

26 06 08 53.6 

B 06 08 52.4 

27 06 08 53.6 

28 06 08 55.4 

31 06 09 03.9 

32 06 09 06.4 

34 06 09 20.8 

35 06 09 30.7 

40 06 10 51.9 

B 06 10 51.8 

42 06 12 34.1 

45 06 12 53.3 

47 06 12 54.1 

B 06 12 55.9 

48 06 13 11.0 

B 06 13 08.8 

49 06 13 29.3 

B 06 13 26.0 

C 06 13 32.2 

50 06 13 47.1 

B 06 13 45.4 

C 06 13 49.6 

51 06 13 59.3 

52 06 14 09.7 

54 06 14 24.9 

B 06 14 24.1 

C 06 14 23.7 



20 39 46.1 
20 39 31.5 

20 36 36.1 

21 22 10.4 
20 36 32.9 
20 42 03.3 
20 34 38.7 
20 26 44.5 
20 39 05.6 
20 40 28.5 

20 37 56.9 

21 31 03.1 
20 36 19.7 
20 33 59.7 
20 32 29.3 
20 37 53.1 

20 38 36.4 

21 31 49.5 
21 33 41.4 
21 38 33.8 
21 37 57.2 
20 29 25.5 
20 40 50.4 

20 42 14.7 

21 50 41.0 

20 39 00.6 

21 23 35.3 
20 29 33.7 
20 30 34.5 

17 54 43.3 

18 00 14.3 
17 59 31.1 
17 57 42.4 
17 58 46.6 
17 58 21.6 
17 55 33.2 
17 56 08.5 
17 56 11.7 
17 54 55.3 
17 55 10.9 
17 54 29.4 
17 52 52.8 
17 43 50.3 
17 45 03.9 
17 44 32.9 
17 43 57.3 



8.60 (0.05) 

7.43 (0.02) 
10.17 (0.02) 

2.48 (0.02) 
10.12 (0.02) 
8.69 (0.01) 
8.18 (0.05) 

5.49 (0.01) 
9.11 (0.01) 
7.82 (0.02) 
7.48 (0.01) 

2.60 (0.01) 
7.69 (0.01) 

8.44 (0.02) 
8.36 (0.02) 

8.61 (0.01) 
7.51 (0.01) 

1.89 (0.01) 
2.66 (0.01) 
3.03 (0.01) 
3.51 (0.02) 
8.35 (0.01) 

7.82 (0.03) 
9.39 (0.03) 
-0.53 (0.01) 
9.34 (0.04) 
-1.00 (0.05) 

7.83 (0.04) 
8.79 (0.02) 
7.71 (0.02) 

8.66 (0.02) 
7.20 (0.02) 
6.25 (0.01) 
5.86 (0.07) 

6.62 (0.03) 
7.85 (0.03) 
8.03 (0.01) 
8.34 (0.03) 
9.18 (0.01) 
9.22 (0.01) 

9.63 (0.01) 
8.94 (0.01) 
6.65 (0.01) 

7.90 (0.02) 

7.67 (0.01) 
7.67 (0.04) 



0.72 
0.58 
0.49 
0.49 
0.51 
0.48 
1.13 
0.42 
0.87 
0.54 
0.78 
0.92 
0.71 
0.91 
0.76 
1.01 
0.55 
0.58 
0.38 
0.93 
1.30 
0.58 
0.55 
0.39 
1.27 
0.64 
0.44 
0.81 
0.52 
0.50 
1.31 
1.07 
0.45 
1.02 
0.42 
0.54 
0.29 
0.37 
0.44 
0.35 
0.36 
0.36 
0.35 
0.53 
0.61 
0.47 
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(0.03 


0.95 


0.16) 


(0.03 


0.94 


0.15) 


(0.02 


1.50 


0.15) 


(0.03 


0.74 


0.15) 


(0.01 


3.43 


0.16) 


(0.01 


2.73 


0.16) 


(0.01 


2.63 


0.15) 


(0.01 


2.01 


0.15) 


(0.02 


1.52 


0.15) 


(0.02 


1.69 


0.15) 


(0.01 


1.88 


0.15) 


(0.04 


0.70 


0.15) 



2.39 


(0.22) 


< 0.79 


< 0.44 


< 0.86 


< 0.48 










5.06 


(0.21) 


0.59 


(0.15) 


0.59 


(0.08) 


< 0.81 


< 0.45 






3.18 


(0.22) 


< 0.80 


< 0.45 


< 0.83 


< 0.46 






4.89 


(0.22) 


< 0.83 


< 0.46 


< 0.85 


< 0.47 






2.39 


(0.11) 


0.29 


(0.08) 


0.31 


(0.04) 


< 0.84 


< 0.46 


< 0.44 


< 0.24 


5.15 


(0.12) 


0.43 


(0.08) 


0.47 


(0.05) 


< 0.87 


< 0.48 


< 0.44 


< 0.24 


4.43 


(0.22) 


0.69 


(0.16) 


0.70 


(0.09) 


< 0.80 


< 0.44 










8.38 


(0.22) 


0.69 


(0.15) 


0.72 


(0.08) 


< 0.78 


< 0.43 










16.27 


(0.22) 


1.72 


(0.15) 


4.06 


(0.08) 


1.17 (0.16) 


2.09 (0.09) 










3.78 


(0.21) 


0.66 


(0.15) 


0.83 


(0.08) 


< 0.79 


< 0.44 










18.21 


(0.22) 


1.45 


(0.15) 


2.67 


(0.08) 


0.63 (0.16) 


0.93 (0.09) 










11.35 


(0.22) 


1.55 


(0.16) 


2.82 


(0.09) 


0.72 (0.16) 


1.33 (0.09) 










11.60 


(0.21) 


0.96 


(0.15) 


1.52 


(0.08) 


< 0.77 


< 0.43 










3.52 


(0.11) 


0.42 


(0.08) 


0.57 


(0.04) 


< 0.84 


< 0.46 


< 0.43 


< 0.24 


3.81 


(0.12) 


0.38 


(0.08) 


0.57 


(0.05) 


< 0.82 


< 0.45 


0.37 (0.09) 


0.25 (0.05) 


12.30 


(0.12) 


0.83 


(0.08) 


1.56 


(0.05) 


< 0.87 


< 0.48 


< 0.43 


< 0.24 


14.04 


(0.23) 


1.29 


(0.17) 


1.85 


(0.09) 


0.71 (0.17) 


0.54 (0.09) 






12.13 


(0.23) 


1.56 


(0.16) 


2.21 


(0.09) 


0.64 (0.17) 


0.71 (0.09) 






4.25 


(0.11) 


0.51 


(0.08) 


0.46 


(0.04) 


< 0.87 


< 0.48 


0.32 (0.08) 


0.24 (0.05) 


14.09 


(0.12) 


1.63 


(0.08) 


3.65 


(0.05) 


1.30 (0.17) 


2.37 (0.10) 


0.36 (0.09) 


0.45 (0.05) 


11.36 


(0.23) 


1.13 


(0.16) 


2.84 


(0.09) 


1.02 (0.17) 


2.04 (0.10) 










12.88 


(0.22) 


1.04 


(0.15) 


1.56 


(0.08) 


< 0.79 


< 0.44 










2.82 


(0.22) 


< 0.76 


< 0.42 


< 0.79 


< 0.43 










0.88 


(0.09) 


< 0.32 


< 0.18 


< 0.49 


< 0.27 


< 0.43 


< 0.24 


8.47 


(0.12) 


0.93 


(0.08) 


2.24 


(0.05) 


0.82 (0.16) 


1.54 (0.09) 


< 0.43 


< 0.24 


2.08 


(0.20) 


0.64 


(0.15) 


0.62 


(0.08) 


< 0.77 


< 0.43 










0.59 


(0.12) 


< 0.44 


< 0.25 


< 0.47 


< 0.26 










3.61 


(0.22) 


0.61 


(0.16) 


0.46 


(0.09) 


< 0.85 


< 0.47 










2.87 


(0.22) 


< 0.79 


< 0.44 


< 0.82 


< 0.46 










3.04 


(0.21) 


0.70 


(0.15) 


0.63 


(0.08) 


< 0.79 


< 0.44 










13.90 


(0.22) 


1.46 


(0.15) 


3.94 


(0.09) 


1.23 (0.16) 


2.62 (0.09) 










8.87 


(0.22) 


1.34 


(0.16) 


2.83 


(0.09) 


1.04 (0.16) 


2.29 (0.09) 










6.83 


(0.22) 


0.94 


(0.16) 


0.97 


(0.09) 


< 0.83 


< 0.47 










2.26 


(0.21) 


< 0.76 


< 0.42 


< 0.81 


< 0.45 










2.34 


(0.22) 


< 0.78 


< 0.43 


< 0.82 


< 0.45 










2.56 


(0.21) 


0.70 


(0.15) 


0.50 


(0.08) 


< 0.80 


< 0.44 










2.41 


(0.21) 


< 0.74 


< 0.41 


< 0.77 


< 0.43 










1.63 


(0.20) 


< 0.76 


< 0.42 


< 0.75 


< 0.42 










11.20 


(0.22) 


1.02 


(0.16) 


1.08 


(0.09) 


< 0.80 


< 0.44 










7.20 


(0.22) 


0.85 


(0.15) 


0.83 


(0.08) 


< 0.76 


< 0.42 










7.28 


(0.20) 


0.75 


(0.14) 


0.50 


(0.08) 


< 0.76 


< 0.42 










5.33 


(0.20) 


< 0.73 


< 0.41 


< 0.73 


< 0.41 










3.81 


(0.21) 


< 0.76 


< 0.42 


< 0.79 


< 0.44 










5.80 


(0.20) 


0.88 


(0.15) 


0.83 


(0.08) 


< 0.75 


< 0.41 










6.85 


(0.21) 


1.06 


(0.14) 


1.46 


(0.08) 


0.65 (0.15) 


0.48 (0.08) 










1.10 


(0.21) 


< 0.75 


< 0.41 


< 0.78 


< 0.43 











Note. — Errors arc given in parentheses. 

B and C denote multiple ammonia pointings that fall within a single 1.1 mm source. 

Upper-limits arc Tmb < 5RMS and W < 5RMSAi7/ViV, where Av is the width of a single channel in velocity, and N is the number of pixels over which the average RMS was calculated. 
denotes sources that were not observed in the NH3(4,4) transition. 



to 

CO 
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TABLE 6 
Derived Gas Properties 



ID 
Number'' 



r(l,l) 



^kin 

(K) 



Tex 

(K) 



a 

(km s^-"-) 



crjvT 
(km s^-"-) 



H2O 



CCS 



emission ; 



01 1.56 

03 0.73 

04 0.96 

05 1.40 

06 1.26 

07 1.43 

11 0.40 

12 3.08 

13 1.68 

B 0.85 

16 1.19 

17 0.47 

18 1.20 

19 0.89 

20 0.53 

21 1.16 

B 1.51 

22 1.29 

24 0.54 

26 0.50 

B 0.65 

27 1.95 

28 0.76 

31 1.64 

32 0.19 

34 0.71 

35 1.81 

40 0.64 

B 0.56 

42 0.47 

45 0.70 

47 0.08 

B 1.02 

48 0.05 

B 0.98 

49 1.22 

B 0.22 

C 0.49 

50 2.43 

B 2.06 

C 2.65 

51 2.60 

52 1.82 

54 0.80 

B 1.25 

C 1.77 



(0.83) 
(0.35) 
(0.54) 
(0.38) 
(0.41) 
(0.19) 
(0.48) 
(0.26) 
(0.14) 
(0.50) 
(0.12) 
(0.19) 
(0.18) 
(0.30) 
(0.27) 
(0.11) 
(0.16) 
(0.18) 
(0.21) 
(0.09) 
(0.22) 
(0.17) 
(0.60) 
(0.85) 
(0.14) 
(0.76) 
(1.33) 
(0.53) 
(0.49) 
(0.49) 
(0.17) 
(0.00) 
(0.26) 
(0.01) 
(0.80) 
(0.66) 
(0.54) 
(0.88) 
(0.18) 
(0.25) 
(0.25) 
(0.35) 
(0.46) 
(0.34) 
(0.27) 
(0.94) 



< 16.76 (1.42) 
17.15 (0.65) 

< 13.93 (1.00) 

< 12.49 (0.72) 
17.45 (0.72) 

14.47 (0.32) 
20.45 (0.93) 
13.39 (0.35) 
25.55 (0.22) 
22.34 (0.86) 
18.93 (0.21) 
25.49 (0.33) 
17.99 (0.32) 

19.48 (0.60) 
19.26 (0.55) 

17.54 (0.19) 
17.04 (0.26) 
20.83 (0.30) 
16.44 (0.37) 
27.47 (0.19) 
28.57 (0.37) 
16.04 (0.25) 

< 15.66 (1.09) 

< 19.15 (1.43) 

29.26 (0.32) 
23.78 (1.36) 

< 18.19 (2.19) 
17.65 (1.02) 

< 15.40 (0.92) 
19.80 (0.90) 

29.55 (0.27) 

32.27 (0.34) 
19.02 (0.46) 

< 18.69 (1.43) 

< 18.44 (1.37) 
19.47 (1.15) 

< 15.72 (0.92) 

< 23.30 (1.48) 
14.21 (0.26) 

15.49 (0.36) 

13.11 (0.34) 

< 14.84 (0.48) 

< 15.54 (0.70) 

20.12 (0.60) 
23.20 (0.44) 

< 22.11 (1.48) 



3.50 
6.25 
4.87 
5.01 
3.88 
5.17 
5.81 
5.05 
6.53 
5.02 
8.85 
10.48 
6.83 
4.17 
5.58 
6.14 
7.85 
7.58 
8.41 
11.61 
7.00 
6.54 
4.86 
3.16 
13.17 
4.22 
3.04 
4.96 
6.35 
6.75 
7.72 
32.27 
7.10 
18.69 
4.30 
4.03 
14.18 
5.35 
6.45 
6.06 
5.46 
4.69 
4.62 
6.16 
5.60 
3.47 



(0.34) 
(1.40) 
(0.96) 
(0.50) 
(0.28) 
(0.24) 
(3.28) 
(0.16) 
(0.23) 
(1.08) 
(0.47) 
(2.73) 
(0.47) 
(0.38) 
(1.24) 
(0.24) 
(0.38) 
(0.50) 
(1.89) 
(1.38) 
(1.25) 
(0.22) 
(1.36) 
(0.23) 
(7.34) 
(1.30) 
(0.34) 
(1.52) 
(2.73) 
(3.65) 
(1.00) 
(0.34) 
(0.89) 
(1.43) 
(0.99) 
(0.52) 
(9.50) 
(4.20) 
(0.19) 
(0.28) 
(0.19) 
(0.20) 
(0.34) 
(1.18) 
(0.46) 
(0.30) 



0.24 
0.24 
0.22 
0.21 
0.25 
0.22 
0.27 
0.22 
0.30 
0.28 
0.26 
0.30 
0.25 
0.26 
0.26 
0.25 
0.24 
0.27 
0.24 
0.31 
0.32 
0.24 
0.23 
0.26 
0.32 
0.29 
0.25 
0.25 
0.23 
0.26 
0.32 
0.34 
0.26 
0.26 
0.25 
0.26 
0.23 
0.29 
0.22 
0.23 
0.21 
0.23 
0.23 
0.26 
0.28 
0.28 



(0.005) 
(0.002) 
(0.003) 
(0.003) 
(0.002) 
(0.001) 
(0.003) 
(0.001) 
(0.001) 
(0.003) 
(0.001) 
(0.001) 
(0.001) 
(0.002) 
(0.002) 
(0.001) 
(0.001) 
(0.001) 
(0.001) 
(0.001) 
(0.001) 
(0.001) 
(0.004) 
(0.005) 
(0.001) 
(0.005) 
(0.008) 
(0.004) 
(0.003) 
(0.003) 
(0.001) 
(0.001) 
(0.002) 
(0.005) 
(0.005) 
(0.004) 
(0.003) 
(0.005) 
(0.001) 
(0.001) 
(0.001) 
(0.002) 
(0.002) 
(0.002) 
(0.002) 
(0.005) 



0.71 
0.58 
0.49 
0.48 
0.50 
0.47 
1.12 
0.41 
0.86 
0.53 
0.77 
0.91 
0.70 
0.90 
0.75 
1.00 
0.54 
0.57 
0.37 
0.92 
1.30 
0.57 
0.54 
0.38 
1.26 
0.63 
0.43 
0.80 
0.51 
0.49 
1.30 
1.06 
0.44 
1.02 
0.41 
0.53 
0.28 
0.36 
0.43 
0.34 
0.35 
0.34 
0.34 
0.52 
0.60 
0.46 



(0.11 
(0.04 
(0.06 
(0.04 
(0.04 
(0.02 
(0.11 
(0.02 
(0.02 
(0.05 
(0.02 
(0.03 
(0.03 
(0.05 
(0.04 
(0.02 
(0.02 
(0.02 
(0.02 
(0.02 
(0.05 
(0.02 
(0.07 
(0.07 
(0.03 
(0.09 
(0.11 
(0.09 
(0.06 
(0.05 
(0.04 
(0.04 
(0.02 
(0.14 
(0.07 
(0.07 
(0.04 
(0.07 
(0.01 
(0.02 
(0.02 
(0.02 
(0.03 
(0.04 
(0.03 
(0.08 



N 

N 
N 
N 
N 
N 
N 
N 
Y 
N 
Y 
N 
N 
N 
N 
N 
N 
N 
Y 
Y 
Y 
N 
N 
N 
Y 
N 
N 
N 
N 
N 
Y 
Y 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 



N 
N 



N 
N 
N 



N 
N 



N 
N 



Note. — Errors arc given in parcntlicscs. 

B and C denote multiple ammonia pointings that fall within a single 1.1 mm source. 
Upper-limits to Tj^^j^ are due to a < 5(T detection in the NH3(2,2) transition. 

Either CCS or H2O was observed, but not both. ■ - - marks sources not observed in the given transition. Y denotes detected, N denotes not 
detected. 
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TABLE 7 
Derived Masses and Densities 



ID 


Miso(120") 


Migo(int) 


Mvii.(int) 


np(int) 


Ilea: 


S(int) 


Nnh3 


NH^iint) 


-^NHs 


Number 


(Mq) 


(M©) 


(M©) 


(103cm-3) 


(103cm-3) 


(10-2gcm-2) 


(1014 cm-2) 


(1021 cm-2) 


(10-8) 


01 


69(7) 


73(6) 


150(28) 


2.39(0.33) 


1.32(0.25) 


1.95(0.08) 


0.53(0.14) 


4.17(0.48) 


1.26(0.37) 


03 


108(14) 


121(13) 


96(10) 


4.31(0.64) 


8.29(0.99) 


3.42(0.18) 


2.31(0.31) 


7.32(0.94) 


3.16(0.59) 


04 


134(12) 


115(9) 


57(8) 


7.18(0.96) 


4.82(0.79) 


4.73(0.15) 


0.78(0.09) 


10.11(1.10) 


0.77(0.12) 


05 


137(13) 


126(19) 


67(7) 


4.52(0.63) 


5.35(0.50) 


3.58(0.21) 


4.01(0.63) 


7.65(0.89) 


5.24(0.35) 


06 


113(15) 


71(11) 


49(5) 


8.54(1.54) 


2.07(0.22) 


4.52(0.46) 


2.28(0.41) 


9.66(1.58) 


2.36(0.57) 


07 


124(21) 


106(18) 


60(4) 


4.86(0.97) 


5.29(0.24) 


3.55(0.42) 


3.55(0.32) 


7.59(1.40) 


4.68(0.96) 


11 


139(13) 


99(12) 


436(56) 


1.98(0.32) 


6.65(1.28) 


1.91(0.13) 


0.95(0.03) 


4.08(0.57) 


2.34(0.34) 


12 


172(23) 


197(24) 


59(4) 


4.23(0.68) 


4.35(0.14) 


3.97(0.26) 


7.04(0.38) 


8.49(1.20) 


8.29(1.25) 


13 


373(40) 


465(39) 


304(18) 


5.58(0.74) 


6.80(0.21) 


6.36(0.12) 


7.58(0.33) 


13.59(1.47) 


5.57(0.65) 


16 


502(53) 


292(28) 


153(8) 


14.29(2.02) 


16.32(0.53) 


10.19(0.35) 


6.65(0.37) 


21.79(2.58) 


3.05(0.40) 


17 


463(30) 


468(39) 


244(17) 


15.69(2.46) 


20.09(1.46) 


12.69(0.24) 


3.53(1.08) 


27.13(2.01) 


1.30(0.27) 


18 


386(38) 


429(37) 


208(15) 


4.94(0.66) 


9.38(0.46) 


5.70(0.12) 


4.93(0.45) 


12.19(1.33) 


4.04(0.58) 


19 


143(15) 


143(15) 


291(27) 


2.64(0.39) 


2.68(0.25) 


2.61(0.12) 


2.92(0.46) 


5.58(0.70) 


5.24(1.05) 


20 


149(16) 


126(14) 


187(17) 


2.95(0.45) 


6.10(0.62) 


2.70(0.14) 


0.78(0.44) 


5.76(0.77) 


1.35(0.78) 


21 


560(57) 


779(66) 


560(31) 


3.76(0.50) 


7.46(0.22) 


5.80(0.12) 


6.22(0.35) 


12.41(1.35) 


5.01(0.61) 


22 


608(61) 


365(35) 


109(7) 


8.56(1.21) 


10.67(0.49) 


7.80(0.26) 


4.54(0.40) 


16.68(1.97) 


2.72(0.40) 


24 


204(19) 


125(18) 


32(2) 


8.48(1.13) 


17.81(1.38) 


5.43(0.11) 


1.51(0.37) 


11.62(1.13) 


1.30(0.13) 


26 


492(31) 


556(20) 


277(8) 


13.69(0.72) 


23.37(0.79) 


12.26(0.19) 


4.21(0.52) 


26.22(1.17) 


1.60(0.21) 


27 


471(47) 


853(65) 


194(13) 


3.58(0.46) 


8.44(0.26) 


5.79(0.06) 


6.72(0.31) 


12.38(1.27) 


5.43(0.61) 


28 


117(10) 


150(12) 


116(19) 


2.17(0.28) 


4.63(0.92) 


2.32(0.05) 


0.65(0.06) 


4.97(0.51) 


1.30(0.18) 


31 


56(7) 


58(7) 


42(8) 


2.29(0.37) 


0.69(0.13) 


1.75(0.13) 


0.16(0.06) 


3.75(0.54) 


0.41(0.18) 


32 


149(15) 


164(16) 


459(28) 


5.72(0.99) 


29.86(1.92) 


4.56(0.31) 


1.82(0.03) 


9.76(0.75) 


1.86(0.15) 


34 


176(9) 


319(10) 


231(42) 


1.37(0.15) 


2.69(0.70) 


2.20(0.07) 


0.42(0.02) 


4.70(0.36) 


0.89(0.08) 


35 


30(7) 


28(4) 


22(6) 


14.41(2.66) 


0.49(0.13) 


4.69(0.66) 


0.14(0.02) 


10.04(1.68) 


0.14(0.03) 


40 


105(10) 


143(9) 


186(26) 


4.98(0.61) 


4.67(0.88) 


3.98(0.06) 


0.80(0.05) 


8.50(0.82) 


0.94(0.11) 


42 


46(11) 


32(8) 


44(6) 


4.81(1.32) 


9.35(1.77) 


2.37(0.52) 


0.65(0.03) 


5.08(1.34) 


1.28(0.34) 


45 


509(53) 


397(32) 


456(28) 


16.96(2.21) 


10.00(0.57) 


12.65(0.22) 


5.74(0.81) 


27.05(2.86) 


2.12(0.38) 


47 


481(49) 


540(43) 


552(37) 


3.90(0.51) 




5.26(0.08) 


2.23(0.03) 


11.25(1.18) 


1.98(0.21) 


48 


84(9) 


87(8) 


310(53) 


2.67(0.38) 




2.23(0.09) 


0.68(0.05) 


4.76(0.57) 


1.43(0.20) 


49 


138(11) 


214(13) 


126(21) 


2.04(0.24) 


2.30(0.41) 


2.51(0.00) 


0.51(0.14) 


5.36(0.50) 


0.95(0.28) 


50 


191(28) 


233(32) 


64(4) 


4.74(0.80) 


8.54(0.24) 


4.53(0.35) 


6.78(0.28) 


9.69(1.47) 


7.00(1.10) 


51 


54(20) 


44(13) 


29(2) 


2.92(0.90) 


3.52(0.18) 


1.89(0.03) 


4.29(0.28) 


4.04(1.20) 


10.60(3.22) 


52 


56(17) 


51(13) 


23(3) 


6.47(1.81) 


3.56(0.33) 


3.36(0.75) 


2.85(0.34) 


7.19(1.93) 


3.96(1.16) 


54 


245(23) 


371(27) 


128(13) 


3.25(0.41) 


7.29(0.81) 


4.11(0.05) 


2.17(0.65) 


8.78(0.88) 


2.47(0.78) 



Note. 



Errors arc given in parentheses. 
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TABLE 8 
Statistical Summary 












units 




Full Sample 






Subset 1 




Property 


minimum 


mean" 


median 


maximum 


minimum 


mean" 


median 


maximum 


J^maj 


pc 


0.18 


0.35(0.10) 


0.34 


0.55 


0.22 


0.37(0.10) 


0.35 


0.55 


^Tfiin 


pc 


0.13 


0.23(0.05) 


0.23 


0.37 


0.16 


0.24(0.05) 


0.23 


0.37 


Rob] 


pc 


0.20 


0.56(0.19) 


0.51 


0.99 


0.30 


0.60(0.19) 


0.56 


0.99 


5^(120") 


Jy 


0.46 


4.30(4.34) 


2.39 


15.38 


0.80 


5.49(4.51) 


3.85 


15.38 


5^(int) 


Jy 


0.43 


4.56(4.62) 


2.09 


17.26 


0.56 


5.84(4.78) 


4.65 


17.26 


7;nb(l,l) 


K 


0.33 


1.76(0.96) 


1.69 


4.15 


0.80 


2.00(0.97) 


1.75 


4.15 


T^b(2,2) 


K 


0.27 


0.82(0.41) 


0.69 


1.72 


0.29 


0.92(0.44) 


0.83 


1.72 


Tex(l,l) 


K 


3.04 


7.27(5.42) 


6.14 


32.27 


3.88 


7.78(5.60) 


6.45 


32.27 


Tk 


K 


12.49 


19.39(4.89) 


18.69 


32.27 


13.39 


20.56(5.08) 


19.47 


32.27 


(Tk'T,.) 


K 


0.00 


12.13(4.61) 


13.05 


21.83 


0.00 


12.78(4.48) 


13.24 


21.83 


VLSR 


km s^-"^ 


-1.00 


6.83(2.96) 


7.85 


10.17 


-0.53 


6.93(2.76) 


7.85 


10.12 


Uv 


km s-i 


0.35 


0.69(0.27) 


0.58 


1.31 


0.38 


0.75(0.27) 


0.71 


1.31 


a 


km s^-"- 


0.21 


0.26(0.03) 


0.26 


0.34 


0.22 


0.27(0.03) 


0.26 


0.34 


O'NT 


km s^-"- 


0.34 


0.68(0.27) 


0.58 


1.30 


0.37 


0.74(0.27) 


0.70 


1.30 


r(l,l) 




0.05 


1.12(0.72) 


1.16 


3.08 


0.08 


1.02(0.70) 


0.80 


3.08 


Mi,„(120") 


M© 


31 


230(180) 


150 


610 


47 


280(180) 


190 


610 


Miso(int) 


M© 


28 


250(210) 


150 


850 


33 


300(220) 


230 


850 


M„ir(int) 


Mq 


22 


190(160) 


150 


560 


33 


220(160) 


190 


560 


M„ir/Miso{mt) 




0.23 


1.00(0.95) 


0.70 


4.37 


0.23 


0.96(0.93) 


0.65 


4.37 


rip (int) 


10^ cm-^ 


1.37 


5.91(4.29) 


4.74 


16.96 


1.37 


6.23(4.42) 


4.81 


16.96 


S(int) 


gcm~2 


0.018 


0.048(0.030) 


0.041 


0.127 


0.019 


0.055(0.032) 


0.045 


0.127 


riex 


10^ cm~3 


0.49 


7.96(6.79) 


6.65 


29.90 


2.07 


9.60(7.05) 


8.29 


29.90 


np(int)/nex 




0.03 


1.50(1.00) 


1.25 


5.22 


0.24 


1.73(1.00) 


1.80 


5.22 


Nnh^ 


IQl^ cm-2 


1.41 


29.68(23.61) 


23.13 


75.76 


4.18 


34.74(23.93) 


29.24 


75.76 


AfHjint) 


1021 cm-2 


3.75 


10.27(6.52) 


8.78 


27.13 


4.08 


11.71(6.94) 


9.69 


27.13 


NNHjN„^{mt) 


io-» 


0.14 


3.00(2.41) 


2.34 


10.60 


0.89 


3.08(2.04) 


2.36 


8.29 



Standard deviation is given in parentheses. 



